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Abstract 

We consider low-temperature behavior of weakly interacting electrons in disordered 
conductors in the regime when all single-particle eigenstates are localized by the 
quenched disorder. We prove that in the absence of coupling of the electrons to 
any external bath dc electrical conductivity exactly vanishes as long as the tem- 
peratute T does not exceed some finite value Tc. At the same time, it can be also 
proven that at high enough T the conductivity is finite. These two statements imply 
that the system undergoes a finite temperature Metal to Insulator transition, which 
can be viewed as Anderson-like localization of many-body wave functions in the 
Fock space. Metallic and insulating states are not different from each other by any 
spatial or discrete symmetries. 

We formulate the effective Hamiltonian description of the system at low energies 
(of the order of the level spacing in the single-particle localization volume). In the 
metallic phase quantum Boltzmann equation is valid, allowing to find the kinetic 
coefficients. In the insulating phase, T < Tc, we use Feynmann diagram technique 
to determine the probability distribution function for quantum-mechanical transi- 
tion rates. The probability of an escape rate from a given quantum state to be finite 
turns out to vanish in every order of the perturbation theory in electron-electron in- 
teraction. Thus, electron-electron interaction alone is unable to cause the relaxation 
and establish the thermal equilibrium. 

As soon as some weak coupling to a bath is turned on, conductivity becomes finite 
even in the insulating phase. Moreover, in the vicinity of the transition temperature 
it is much larger than phonon-induced hopping conductivity of non-interacting elec- 
trons. The reason for this enhancement is that the stability of the insulating state 
is gradually decreasing as the transition point is approached. As a result, a single 
phonon can cause a whole cascade of electronic hops. 
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1 Introduction 

Transport properties of conducting materials at low temperature T are deter- 
mined by an interplay between the interaction of the itinerant electrons with 
each other and the quenched disorder which creates a random potential act- 
ing on these electrons. In the absence of the electron-electron interaction all 
physics is dominated by the phenomenon of Anderson localization [1] - e.g., dc 
electrical conductivity a can be qualitatively different depending on whether 
one-particle wave functions of the electrons are localized or not. In the latter 
case cr(T) has a finite zero-temperature limit, while in the former case o-(T) 
vanishes when T ^ 0. Therefore, Anderson localization of electronic states 
leads to the Metal to Insulator Transition at zero temperature. 

When speaking about zero temperature, we need to consider only localization 
of the electronic states close to their Fermi level. The conductivity becomes fi- 
nite at any finite temperature provided that extended states exist somewhere 
above the Fermi level. It is commonly accepted now that localized and ex- 
tended states in a random potential can not be mixed in the one-electron spec- 
trum and thus this spectrum in a very general case is a combination of bands 
of extended states and bands of localized states. A border between a localized 
and an extended band is called mobility edge. If the Fermi level is located 
inside a localized band and inelastic scattering of the electrons are completely 
absent, the conductivity should follow Arrhenius law cf{T) oc ex.p{—Ec/T), 
where E^. is the distance from the Fermi level to the closest mobility edge. 
Another common belief following from the scaling theory of Anderson local- 
ization [2,3], is that in low dimensionality d, namely at d = 1,2 all states are 
localized in an arbitraryly small disorder, while for free electrons (no periodic 
potential) > is finite at d — 3 (for d — 1 this statement was proven 
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rigorously both for thin [4,5] and thick [6,7] disordered wires). It means that 
without inelastic processes 0-^=1,2 (7") = 0, while for (7^=3 (T) one should expect 
the Arrhenius law. Note that for electrons in a crystal within a given con- 
duction band the latter conclusion could become incorrect - strong enough 
disorder can locahze the whole band. 

Situation becomes more complicated when the inelastic processes are taken 
into account. In particular, electron-phonon interaction leads to the mecha- 
nism of conductivity known as hopping conductivity [8,9] - with an assistance 
of phonons electrons hop between the localized states without being activated 
above the mobility edge. As a result, o'(T) turns out to be finite (although 
it can be very small) at any finite T even when all one-electron states are 
localized. 

Can interaction between electrons play the same role and cause the hopping 
conductivity? This question was discussed in literature for a long time [10- 
14] and no definite conclusion was achieved. The problem is that although 
the electric noise exists inside the material with a finite ac conductivity ^ the 
"photons" in contrast with phonons become localized together with electrons. 

In this paper we demonstrate that electron-electron interaction alone cannot 
cause finite conductivity even when temperature is finite, but small enough. In 
the absence of phonons and extended one-electron states a system of interact- 
ing electrons has exactly zero conductivity below some temperature T^.. This 
critical temperature is infinite if the distance between the electrons exceeds 
the localization lengths of all electronic states. In the opposite case Tc turns 
out to be finite and depends on the typical number of electrons within the 
localization volume as well as on the strengths of the electron-electron inter- 
action. We also argue that at high temperatures T > T^. the conductivity cr{T) 
is finite. It means that at T = Tc the system of interacting electrons subject to 
a random potential undergoes a genuine phase transition that manifests itself 
by the emerging of a finite conductivity! 

This transition can be thought of as many-body localization - it applies to 
many-body eigenstates of the whole system. This localization occurs not in the 
real space, but rather in the Fock space. This fact does not affect the validity 
of the concept of mobihty edge. In fact, the existence of the "metaUic" state 
at T > Tc implies that the many-body states with energies £ above £c are 
extended. One can estimate the difference between Sc and the energy of the 
many-body ground state Sq as Sc — Sq ~ TM{T)^ where M{T) is the total 
number of one-particle states in the energy strip of the width T. Note that 
the existence of the extended many-body states above the mobility edge does 
not contradict the fact that below Tc there is no conductivity - in contrast 

^ In this paper we mostly focus on dc conductivity. As to ac conductivity, it never 
vanishes, because at any frequency density of resonant pairs of states is finite. 
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with the case of onc-particle locahzation there is no Arrhenius regime since 
— Sq turns out to be proportional to the volume of the system, i.e., is 
macroscopically large (see Sec. 2.1.2 below for more details). 

In order to avoid possible misunderstanding we would like to emphasize at 
the very beginning that throughout this paper we will be focused only on the 
inelastic collisions between the electrons, i.e.. on creation or annihilation of 
real electron-hole pairs. There are other effects of electron-electron interac- 
tions - they can be understood as renormalization of the one-particle random 
potential by the interaction. This renormalization is temperature dependent 
and thus leads to a number of interesting effects, for example, the interaction 
corrections to the density of states and conductivity in disordered metals [15]. 
On the insulating side of the one-particle localization transition similar effects 
cause the well-known Coulomb gap [9] which suppresses hopping conductivity. 
On the other hand, this is still a time-independent correction to the time- 
independent random potential. As such, it can maybe shift the position of the 
many-body Metal to Insulator transition, i.e., renormalize Tc, but is unable 
to destabilize the insulating or metallic phases. From now on we will simply 
neglect all elastic (Hartree-Fock) effects and concentrate on the real inelastic 
electron-electron coUisions. 

Localization of the many-body states in the Fock space has been discussed in 
Ref. [16] for the case of zero-dimensional systems with finite, although maybe 
large, number of electrons. Authors of Ref. [16] used an approximate mapping 
of the Hamiltonian of a metallic grain with large Thouless conductance g and 
moderate interaction between the electrons to the one-particle Hamiltonian 
on a lattice with the topology of the Cayley tree and an on-site disorder. 
The latter problem has an exact solution [17] that contains the localization 
transition. In terms of interaction electrons this transition means that below 
certain energy of the excitation the one-particle excitation states are quite 
close to some exact many-body excitations. As to the one-particle excitations 
with energies higher than the critical one, its wave function can be viewed as 
a wave packet, which involves a large number of the many-body eigenstates. 
Being a property of a finite system this transition could be nothing but a 
crossover, which becomes more pronounced with increasing of g. 

For an infinite system with d > the presence of spatial degrees of freedom 
makes the situation more complex. In this case, Cayley tree approximation 
turns out to be insufficient. Nevertheless, a consistent analysis of a model with 
weak and short range interaction to all orders of perturbation theory enabled 
us to analyze the many-body localization transition and to demonstrate that 
both the metallic state at high temperatures and the insulating state at low 
temperatures are stable and survive all higher loop corrections to the locator 
expansion. This allows us to claim that the existence of the transition is proven 
on the physical level of rigor. 
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It should be noted that such an insulating state that is characterized by exactly 
zero conductivity is quite different from all other known types of Metal to 
Insulator transitions. For example, Mott insulator is believed to have finite, 
though exponentially small conductivity at finite temperatures. 

The remainder of the paper is organized as follows. Sec. 2 contains the dis- 
cussion of the problem on the qualitative level: we define the many-body lo- 
calization, show its macroscopic implications, and discuss the relation to the 
Anderson model on a certain graph. In Sec. 3 we specify the model many-body 
system to be studied throughout the paper. In Sec. 4 we show how to treat this 
model in the framework of non-equilibrium Keldysh formalism, introduce the 
self-consistent Born approximation (SCBA), and derive the general equations. 
Using these equations, we demonstrate the existence of the metallic state at 
high temperatures and study its properties in Sec. 5. Sec. 6 is devoted to the 
proof of the existence of the insulating phase at low temperatures; we evalu- 
ate the transition temperature as the limit of stability of the insulating phase. 
Sec. 7 is dedicated to justification of SCBA; we demonstrate that corrections 
to SCBA are indeed small. Finally, in Sec. 8 we summarize the results and 
present an outlook of the future developments. 



2 Qualitative discussion 

2.1 Macroscopic manifestations of the many-body localization transition 
2.1.1 Single-particle localization 

Let us begin with the brief review of the basic notion emerged in a study of 
the properties of the one electron wave functions in a disordered potential in 
d dimensions. Depending on the strength of the disorder potential, a wave 
function (pair) of an eigenstate a corresponding to a given energy can be 
classified as localized or extended: 



where Qoc is the localization length which depends on the eigenenergy ^q,, 
and V is the volume of the system. Exponentially locahzed states have a 
maximum amplitude at some point in the system and a wave function 
envelope which falls off exponentially, whereas the extended states spread 
relatively uniformly over the whole volume of the system. It is believed that 
the coexistence of the localized and extended states in the same energy range 
is not possible, and the spectrum splits into bands of localized and extended 




extended. 



localized; 



(1) 
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states. The energies separating such bands are known as mobility edges. For 
example, for free electrons in a disorder potential in the dimensionality three 
and higher, only one mobility edge £i exists, so that 



< '■ localized; 

> £i : extended. (2) 

The statement about the asymptotic behaviour of the single-particle wave- 
functions (1) can be translated into the property of the matrix elements of 
a certain local operator A{R) (it might be the local mass or current density 
operator, etc.): 

A^p{R) = j d''r<j>l{r)A{R)<t>p{r). (3) 

Then, 

{ \r\ 
< e '^loc ■ locahzed; 
^ [it, ) ' ^^t^^d^d' 

where the linear scale L^^ is controlled by the transmitted energy cua/s — ^a — Ci3 
only, and L^^ — oo for — > 0. No summation over the repeating indices is 
implied in Eq. (4). Energies and are assumed to be sufficiently close to 
each other, so that the localization length is approximately the same for the 
two states. 

The notion of the matrix elements (4) is intimately related to the observable 
quantities. As an example, consider the Kubo formula for the density-density 
response function 

where the overlap £^^(r) is given by Eq. (4) built on the operators of the local 
density, Q^R), such that Q{R)\r) = 6{R — r)\r). In the thermal equilibrium 
for the fermionic particles, the occupation numbers are given by the Fermi 
distribution. 

In a conducting system, the low-frequency asymptotic behavior of the density- 
density response function is always determined by diffusion as the total number 
of particle is the only conserved quantity in the system: 

n(u;;r)oc /d'^Q— ^J^^, (6) 
J —luj + DQ 
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where D is the diffusion coefficient. First, let assume that all the states are 
localized. Then, comparing Eq. (5) with Eq. (4), we see that the asymptotic 
behaviour of the wave function precludes the diffusion propagation (6). There- 
fore, the diffusion coefficient vanishes for any temperature T. The same is true 
for the dc conducivity, as it is related to the diffusion coefficient via Einstein 
relation, a cc D. Physical meaning of vanishing D and a is the impossib- 
lity for an excitation caused by a local external perturbation to propagate 
in a localized system, and uniformly span all the phase space allowed by the 
conservation of energy. 

If a finite mobility edge (2) exists and the Fermi level ep lies in the energy 
region of localized states, the conducitivity is determined by the exponentially 
small occupation number of the delocalized states 

a{T) oc e-(^i-^^)/^. (7) 

In this paper we will be interested in transport properties of the systems 
where all single-particle states are localized, and thus without many-body 
effects o" = at any temperature. It is well established now that the mobility 
edge does not exist for one- and two-dimensional systems for any disorder, and 
all single-particle states are indeed localized. Such a situation can also arise 
in an arbitrary high dimensionality as well. Indeed, consider as an example 
Anderson model with one state per lattice site [1]. It is well known that all 
eigenstates become localized as soon as the on-site disorder exceeds a critical 
value. 



2.1.2 Many-body localization 

Let us now turn to the discussion of the many-body localization. From now 
on, we assume that all one-particle state are localized in the sense of the 
previous subsection. Can the interaction cause finite conductivity? Consider a 
many-body cigcnstate of the interacting system, with the corresponding 
eigenenergy Ek- Our purpose is to generalize the notion of localization to such 
many-body states. 

In the coordinate representation, the many-body wave function ({^jjjLi) 
depends on the coordinates of all particles in the system. The single-particle 
states forming this many-body state are located everywhere in the volume V. 
Thus, no definition analogous to Eq. (1) can be constructed. On the other 
hand, the relation (4) can be used, if one takes a local additive one-particle 
operator: 

A(J^) = 5] A,^(i^)4c^ = / d'rc^{r)A{R)c{r), (8) 

ad •' 

where we introduced the fermionic creation and annihilation operators cj^, 
c^(r), Cq, c(r) in the basis of the one-particle eigenstates and in the coordinate 
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representation, respectively. 



Wc consider the matrix elements of the local operator between two exact 
many-body eigenstates: 



Akk'{R) = A{R) ^k') ■ (9) 
Then we can define localized states by a relation, analogous to Eq. (4): 

< e «(^fe) ; locahzed; 
^(4;); extended. 

Again, uO}^j.i = E^ — Ey and the energies i?^ and Ei^i are taken sufficiently close 
to each other, so that the difference between C(-^fc) and Ci^k') need not be 
taken into account. ^ 

Equation (10) is the definition of the localized many-body states. To elucidate 
its meaning, let calculate the matrix element (9) between two Hartree-Fock 
states 

= n (l - /a + faci) Ivacuum) , (11) 

a 

where the state \^k^) is completely characterized by the set of occupation 
numbers = 0, 1, corresponding to empty and filled single-particle states, 
respectively. Substituting Eq. (11) into Eq. (8) we find for k ^ k' 



^kk' — XI 



n (12) 

where we omitted the argument R in both sides of the equation. The states 
k and k' connected by the operator A[R) are obtained from each other by 
creation of a single electron-hole pair, so that only one term in the sum over 
a and is different from zero. We note that the matrix elements Aa/siR) 
are exponentially suppressed unless both states a and [3 are located near the 
point R. Thus, the distance between the electron and the hole cannot exceed 
the single-particle localization length (loc, and the number of states k' which 
can be connected to the given state k by the local perturbation is effectively 
finite, even though the total number of the Slater determinat states (11) scales 
exponentially with the size of the system. 



^ Let us emphasize that the criterion (9) can be extended to the arbitrary number 
of electron-hole excitations, A{R) Ai{R)A2{R) . . .An{R), corresponding to the 
local heating of the system. The asymptotic behaviour will still be determined by 
Eq. (10). 
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We substitute Eq. (12) into the localization criterion (10) and obtain 

Because of the (5-symbols, \Ek - Ek'\ = \ J2^,^^{f!^ - /^')| = - ^/sl, where 
a and f] arc the only two states contributing to the sum for given k, k'. If 
ia,i)3 are not too far from the Fermi level, so that Ooc(Ca) ~ CiodC/s), we 
find that Cf,^, (r) has the same spatial dependence as the kernel (4) for non- 
interacting system. Thus, according to definition (10) any Hartree-Fock state 
(11) is localized with the locahzation length ({Ek) — C,ioc- 

Let us apply the same ideas to the Kubo formula for the many-body density- 
density response function 



n (a;; r) = ^ P^n^ (a;, r) = Y. [Pk - Pw] ^kk' (a;, r) (14a) 

k kk' 

Hfe {u, r) = 5Z [nfefe' - n^'fe (a;, r)] ; (14b) 

k' 

ImHfefe' (u;, r) = ^>C^fc,(r)5 {u + Ek- E^') , (14c) 

where Pk is the probability for the system to be in the eigenstate k and 
Ilfc {uj, r) is the linear response of the system in this particular eigenstate. 
The overlap function £^^,(r) is given by Eq. (10) with the operator of local 
density g{R) = c^{R)c{R). 

Formulas (14) are valid for any many-body eigenstates l^'fe). For the Hartree- 
Fock state (11), we can substitute Eq. (13) into Eqs. (14) and obtain Eq. (5) 
with 

/a = E Pj^fa- (15) 

k 

For the initial Gibbs distribution 

Pk = ^""^ ^^^i N , (16) 
Efcexp(-f) 

and the Hartree-Fock spectrum, Ek — J2a fa^a, formula (15) gives the Fermi 
distribution function. 

According to Eqs. (13) and (14c), the spatial dependence of the correlator 
H(c<j, r) is still given by Eq. (4) for the localized single-particle wave functions. 
Therefore, both the conductivity and the diffusion coefficient vanish for the 
Hartree-Fock states constructed in the basis of exact one-particle wave func- 
tions of disordered systems. The same statement can be made about certain 
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wave functions which are formally very different from a single Slater determi- 
nant. Namely, consider the creation of the electron-hole pair on top of some 
eigenstate 1^^^) and expand the result in terms of other eigenstates 



E 

k' 



r-ikk' 

rap 



1. 



(17a) 



For the Hartree-Fock state (11), only one term contributes to the sum, see 
Eq. (13). It is possible to check that the state would remain localized in a 
sense of Eq. (10), if the number of terms contributing to the sum in Eq. (17a) 
is large but finite, i.e. 



lim 

V-*-oo 



E 

. k' 



\/^kk' 
ral3 



< OO. 



(17b) 



In complete analogy with the non-interacting system, this corresponds to the 
insulating or localized many-body state; excitation can not propagate over all 
states allowed by the energy conservation. 

Conductivity can be different from zero only if the wave functions of the 
excitations can be broken onto the infinite number of eigenstates 



lim 

V^oo 



EI^J 

. k' 



kk' 

:/3 



= OO, (17c) 
which would correspond to the metallic or extended many-body state. 



Developed metallic state corresponds to the case when electron-electron in- 
teraction mixes the excited state with all the eigenstates in the system with 
close enough energy: 



IC^jToc "5(Efe + o;„^-i?fcO 



(17d) 



where 5-function should be understood in the thermodynamic sense: its width, 
although sufficiently large to include many states, vanishes in the limit V 
OO. Only in this regime, which may also be called ergodic many-body state, the 
electron-electron interaction can bring the system from the intitial Hartree- 
Fock state to the equilibrium correspoding to spanning all the states permitted 
by the energy conservation. In this case, the averaging over the exact many- 
body eigenfunction is equivalent to averaging over the microcanonical distri- 
bution, and temperature T can be defined as a usual Lagrange multiplier. It 
is related to Ek by the thermodynamic relation: 



J. 

Ek-Eo^ I Cv{T,) dTi, 



(18) 



where Eq is the ground state energy, and Cy(T) oc V is the specific heat. 
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The main result obtained in the present paper is the proof of existence of 
the extensive many-body mobility edge Sc oc V. This proof is based on two 
statements: (i) states with sufficiently large energies — Eq are extended; 
(ii) states with sufficiently small energies E^ — Eq are localized. 

The first statement follows from the validity of high-temperature expansion 
for the quantum corrections to conductivity [15]. The quantum corrections to 
conductivity are divergent for o? = 1, 2; for non- interacting systems these weak 
localization corrections have no cutoff other than the size of the system. In 
interacting systems two kinds of new phenomena appear: (a) interference due 
to the scattering off the self-consistent Hartree-Fock potential (see Ref. [18]), 
and (b) inelastic electron scattering. Effects of type (a) are regularized by the 
temperature itself and do not produce any consequences for the present paper. 
At the same time, inelastic electron scattering leads to regularization of the 
weak localization corrections due to appearence of the inelastic rate and the 
corresponding length: 



where A < 1 is the dimensionless interaction constant, D is the diffusion 
coefficient, and g{L) is the dimensionless conductance of the li-dimensional 
cube of linear size L. All the interference corrections are finite if 



Using g{Cioc) — 1, we rewrite Eqs. (19) and (20) as 



Inequality (21) is the condition of appUcability of the expansion from the 
metallic state. In fact, it is also valid for higher-dimensional systems with the 
finite bandwidth E^ and all single-particle eigenstates localized, if Ei, ^ T*^^'). 

At T < Tel the perturbation theory breaks down. It may indicate either (i) a 
simple insufficiency of the perturbation theory to describe the metallic state, 
or (ii) existence of the many-body mobility edge Ec- It important to emphasize 
that £c is an extensive quantity £c oc V. 

In case (i), states are extended for all E^, so that the conductivity would 
remain finite (no matter how small) down to zero temperature. In case (ii) the 
many-body eigenstates with Ek — Eq < £c are localized, see Eq. (17b). 

In case (ii), the partial conductivity of one state, 0"^, defined analogously to 
Ilfc from Eqs. (14), is zero for Ek — Eq < Ec- 




(19) 




(20) 
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Fig. 1. Schematic temperature dependence of the dc conductivity cr(T). Below the 
point of the many-body metal-insulator transition, T < Tc, cr{T) = 0, as shown in 
Sec. 6. Temperature interval T > T^"*) > Tc corresponds to the developed metallic 
phase, where Eq. (17d) is valid. In this regime for the model described in Sec. 3 
cr(r) is given analytically by Eqs. (93)-(99) and plotted on Fig. 10. At T > T^^^) 
the high-temperature metallic perturbation theory of Ref. [15] is valid. 

Results of Rcf. [16] strongly indicate that the case (ii) is realized. We will 
review, extend and refine the arguments of Ref. [16] in the next subsection; 
here, we simply proceed with the discussion of the macroscopic manifestations 
of this scenario. Let us assume that the equilibrium occupation is given by the 
Gibbs distribution (16). One could think that it would still imply the Arrhenius 
law (7) for the conductivity. However, this is not the case for the many-body 
mobility threshold. In fact, in the limit V — > oo 



C7(T)=0; T<Tc, 
where the critical temperature is determined by Eq. (18): 



dT,CviT,)=£c. 



(22a) 



(22b) 



The schematic temperature dependence of the conductivity is summarized on 
Fig. 1. Therefore, the temperature dependence of the dissipative coefficient in 
the system shows the singularity typical for a phase transition. 

To prove Eqs. (22) we use the Gibbs distribution and find 

a(T) =}_^P,a{E,) = jo. ^^^s(e)-e/t ' 

where the entropy S{E) is proportional to volume, and E is counted from the 
ground state. The integral is calculated in the saddle point or in the steepest 



13 



decent approximations, exact for V — > oo. The saddle point E{T) is given by 



dS_ 
dE 



E=E{T) 



T 



1 



Taking into account cr{E) = for £■ < we find 




As both energies entering the exponetial are extensive, E{T),Sc oc V, we 
obtain Eqs. (22). 

As we aheady mentioned, vanishing of the dissipative conductivity at T < Tc 
means freezing of all relaxation processes. In particular the micro canonical 
distribution could never be established for the closed system. In this respect, 
the dynamics of the system resembles the glassy state [19]. 

To establish the thermal equilibrium in such insulating state requires finite 
coupling of the system with the external reservoir (i.e., phonons). The presence 
of the finite electron-phonon interaction (as phonons are usually delocalized), 
smears out the transition, and o"(T) > for any temperature. Nevertheless, 
for the weak electron-phonon interaction, the phenomenon of the many-body 
metal-insulator transition (22) manifests itself as a very sharp crossover from 
phonon induced hopping at T < Tc to the conductivity independent of the 
electron-phonon coupling at T > Tc. 



2.2 Microscopic mechanism of the many-body localization transition 

As discussed in the previous subsection, the existence of extended many-body 
states at high energies is an established fact [15]. Here we focus on the existence 
of localized states at low energies, and discuss the correspondence between a 
many-electron interacting system and the Anderson model on a certain graph. 
As we have already mentioned, it is convenient to analyze many-body localiza- 
tion in terms of single-particle excitations cj^l'^k) above a certain eigenstate 
of the interacting system, namely, how these excitations spread over other 
many-bogy eigenstates (consideration of the electron-hole excitation is per- 
formed in a same fashion and does not bring anything qualitatively different) . 
This discussion has a close relation with that of Ref. [16]. 

Let the system initially be in the eigenstate At time t = an extra 

electron is created in the single-particle state a. The many-body Schrodinger 
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equation describing the subsequent time evolution of such state I'^kai't)) is 

'zdt-H]\^ka{t)) = mci\'^k), (23) 

where the right-hand side determines the initial condition at t — 0, and the 
Hamiltonian, written in the basis of the exact single electron wavef unctions, 
is given by 

1 

^ = ciCa + 2 II 4^^^^^^ . (24) 

The interaction part of the Hamiltonian contains only non-diagonal terms, 
a ^ ^,S; P ^ j,5, and assumed to be antisymmetrized, Vapjs — ~Vi3a-yS — 
—Vaps-y- The diagonal matrix elements are already included into the definition 
of the Hartree-Fock spectrum {^q}. 

We make the Fourier transform of Eq. (23): 

(e + Ek - H)\^ka{.e)) ^ ci\^k) (25) 

(here we count the energy e from that of the reference state E^), and solve 
this equation for |^fea(e)) by iterations: 

\^ka{e)) = (ICH^)) + IV'S(e)) + ■■■)■ (26) 

The zeroth order in the electron-electron interaction is just the one-particle 
excitation itself: 

^S) = 4l*.)- (27) 
The first order corresponds to a three-particle excitation above |^fe): 

4L^)=E ^44^/31*.), (28) 

where we introduced a shorthand notation 

5^-^7 + ^5-^/3 (29) 
for the energy of the three-particle excitation. 

As usual in the theory of metal-insulator transition [1], one starts from esti- 
mating the probability for the sum (28) not to be small. We notice that the 
geometric distance between particles is of the order of the one-particle local- 
ization length Qoc-i as the interaction is short-ranged. Thus, the denominators 
in Eq. (28) are the random quantities with the maximal value of the order 
of the level spacing in the localization length 5q. Assuming weak interaction, 
A = max|Va/375|/5(^ <^ 1, we see that the value of sum is, in fact, determined 
by the smallest denominator. 
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Let us assume that the typical number of terms {i.e., terms not suppressed by 
the matrix elements) is 3> 1. As the denominators are distrubited approx- 
imately uniformly for —5,^ < e — S^^ < 5^, the smallest denominator is of the 
order of 5(^/K. Therefore, if K\ <^ 1, the probability to find the large mixing 
is small, whereas for 

K\ > 1 (30) 

it becomes of the order of unity. This estimate up to a factor of the order of 
I In A I is the basis for finding the position of the transition [1]. 

The remaining non-trivial problems are: (i) to find the connectivity K, and 
(ii) to check that higher-order terms indeed match the locator structure of 
Anderson, i.e., the the number Hn of the nth order terms scales as 

Mn ^ K^, (31) 

where the prefactor can be any algebraic function of n. Condition (31) is very 
important: if jVn ^ K"', e.g., Mn — n\, then the probability to find small 
enough denominators in high orders of the perturbation theory Mn^"' — 1 
would be always of the order of unity no matter how small the interaction 
constant A is; the system in this case would always remain metallic. If, op- 
positely, the number of terms in high orders is small, e.g., Mn — K''^ /n\, the 
higher-order terms will not contain the small denominators, even though the 
lowest orders of the perturbation theory do. In this case, the system would 
always remain insulating. In both cases, the conlusions drawn from the lowest 
order perturbation theory [16] would be misleading. 

To estimate K from Eq. (28), one has to account for (i) the structure of the 
eigenfunction \^k)) i-^-i which indices a of fermionic operators are allowed 
not to produce zero result after action on the wave function; (ii) the structure 
of the matrix elements Vap-yS- 

The energy of a typical Hartree-Fock state is almost uniformly distributed 
among the localization volumes. The fact that the eigenstate is a linear 
combination of an infinite number of Hartrcc-Fock states does not affect this 
assumption. Therefore, the state |^fc) is contributed by electron-like and hole- 
like excitations with the energy smaller than T, where T is of the order of 
temeraturc given by Eq. (18). Considering the action of the fermionic operators 
on such a state we estimate 

-T<^^,^8; ^p<T. (32) 

Second restriction is based on the structure of the matrix elements. As we 
already mentioned, they restrict the final states to be at the distance not ex- 
ceeding the localization length C,ioc from each other. In addition to this spatial 
restriction, the matrix elements decrease rapidly when the energy difference. 
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say — exceeds the level spacing 5^, i.e., Vajs-^s — A(5^ only for 

\ia - i5\M^ - Q ^ k or \ia - - is\ < (33) 

and vanish rapidly otherwise. 

Counting the number of states /3, 7, 5 within one localization volume, such 
that they produce the value of the denominator in Eq. (28) smaller than 5^. 
and satisfy both restrictions (32) - (33), we obtain for T ^ 5^ ^ 

K^^. (34) 

k 

Comparing Eq. (34) with the estimate (30) we conclude that at T > ^ the 
first correction to the wavefunction (28) contains at least one term of the order 
of unity, which may signify the transition in the system. 

To check that this conclusion is not an artefact of the lowest order perturbation 
theory, one has to analyze further terms in the expansion (26). The second 
order correction is given by 

c«i/3i /3,7,5 ^ ^ "7<5 

+ 2^ 2^ _^/37i "^ai^^iW^Til^fc/ (35) 

ai,/3i,7i fS,j,S ^ "7ai/3i ^ "7^ 

ai,7i,(Si /3,7,5 ^ "75ai ^ "7<5 

where the (2n + 1) -particle energies are defined as 

n+l n 

^t:t,. = EU-i:^Pr (36) 



i=l i=l 



The first term in Eq. (35) is again a three-particle excitation. We can estimate 
its typical value as A^T, which is smaller than the first-order contribution (28) 
by a factor A <S 1; similar type terms arise also from the second and the third 
lines for c^, c with coinciding indices. The actual meaning of this terms is the 
renormalization of the value of the two-particle interaction strength. They are 
analyzed with more rigor in Sec. 7.1, here we simply neglect them as they do 
not affect the statement about the transition. 



Once again, we estimate the number of the relevant terms in the multiple 
sum in the last two lines Eq. (35). They are shown pictorially on Fig. 2b. 

^ The estimate of Ref. [20] for the quantum dot at T = is ~ e^/^^j where 5 is 
the level spacing in the quantm dot. Our estimate is different because the energy 
restriction on the matrix elements (33) is absent in a quantum dot. 
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(2n+l)! (?) 

Fig. 2. Pictorial representation of the paths corresponding to the first (a), second 
(b), third (c), and higher (d) orders of the perturbation theory for the mixing of 
one-particle excitations with many-particle states. Process (d.2) is obtained from 
the process (d.l) by the application of the interaction operators once to each of 
(2n+l) final particles of the process (d.l). The sequence of such applications can 
be ordered in (2n + 1)! ways, each producing formally different path. 

The structure of the state gives the restriction similar to Eq. (32) on 
the electron (entering Eq. (36) with plus sign) and the hole (entering Eq. (36) 
with minus sign) energies respectively. The matrix element Va/s-ys restricts the 
energies by Eq. (33). Combining those two restrictions and requiring each 
denominator in Eq. (35) to be smaller than the level spacing 5^, one estimates 

M K\ (37) 

where K is given by Eq. (30). Aparently, Eq. (37) agrees with the require- 
ment (31). However some questions may arise already on this stage. 

One notices from Eq. (35) that there are terms which can be obtained from 
each other by the simple permutation of indices of the creation or annihilation 
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operators. The states obtained by such permutations from each other are in 
fact identical, so one may be tempted to amend Eq. (37) as [20] 

n\{n + Ij! 

where the combinatorial factors describe the indistinguishability among n + 1 
electrons and n holes [n — 2 for Eq. (35)]. 

However, in the theory of Anderson transition [1], the relevant parameter 
is not the number of available final states, but the number of statistically 
independent paths leading to these states from the initial state {e.g., for the 
Anderson transition in three dimensions the number of states grows as 
whereas the number of paths grows exponentially with n). By construction, 
the summation in Eq. (35) is performed over the single-particle indices and 
not over the final states. Moreover, the denominators in each terms involve 
different combinations of levels and thus produce statisitically independent 
contributions. This corresponds to the summation over paths number of which 
cancels the factorials in Eq. (38). 

Having demonstrated the absence of the factorial factors suppressing the num- 
ber of terms Hm one may worry about the opposite extreme: factorial growth 
of Mn- The way to obtain this factorial is shown on Fig. 2d. The transition be- 
tween the two states with (2n+ 1) excitations and (6n + 3) excitations may be 
realized in (2n-|- 1)! ways. These ways are different by the order in which inter- 
action operator acts on the particles, thereby producing different intermediate 
states. If these paths were statistically independent, it would mean 

jV3n+l~(2n+l)!K2n+l_^^^ (39^ 

so that no transition would occur and the system would always be delocalized. 

However, such paths are, in fact, correlated and the sum of the correspond- 
ing amplitudes produces always the result of the order of an amplitude of a 
single path; the terms essentially cancel each other. To demonstrate this can- 
cellation we consider the third order of the perturbation theory. The terms 
schematically shown on Fig. 2c correspond to 



/(3)\ _ 2V^^^^^^j / Vg 1/317^1 Kti/3i75i 

X ^2 ^2 ^Li 4i ^/^^-Ji C72 1 ■ 

^ "75 



(40) 



Two terms in the second factor correspond to two ways to obtain the final 
seven-particle state. One immediatley notices that the matrix elements in the 
two terms are the same. Moreover, the denominators in two terms can be 
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combined as 

I _ "^7^ ~ ^ "ai/3ia2/32 ( a-i \ 

Generally, one seeks to maximize the transition amplitude by choosing the 
smallest possible denominators for each of the three factors in Eq. (40) inde- 
pendently. Each denominator being of the order of S^/K, the whole expression 
would be proportional to K^. However, one can see immediately from Eq. (41) 
that small denominators of the first and the third factors of Eq. (40) appear 
in the numerator of the second factor. Thus the amplitudes shown Fig. 2 
(c.2)-(c.3) are not indepenedent, and in fact, cancel each other, producing a 
contribution proportional to instead of K^. One can follow similar cancel- 
lations in any order of the perturbation theory. Therefore, the estimate (39) 
is not valid and scaling (31) remains intact. Therefore, one can use the An- 
derson result [1] for the critical connecitivity and obtain the estimate for the 
transition temperature 

the many-body mobility edge £c is then found from Eq. (22b). The qualitative 
estimate (42) is in agreement with the result of a quantitative calculation 
performed in subsequent sections for a specific model [Eq. (160)]. 

The arguments given in this section show that, even though the many-body 
problem (24) exhibits a close analogy with the problem of Anderson localiza- 
tion [1], an exact mapping to the Anderson model on some graph, like one 
suggested in Ref. [16], is problematic. First, the geometric structure of this 
graph is unknown. Analogy with the Cayley tree, popular due to the exact 
solvabihty of the corresponding Anderson model [17], is not apphcable, strictly 
speaking, because in the many-body problem two states can be connected in 
more than one way, which would be prohibited for the Cayley tree. Second, 
the many-body problem does not allow for simple counting of statistically 
independent paths, which was the main idea of the Anderson's solution [1], 
as transition amplitudes, corresponding to seemingly different paths, exhibit 
striking correlations. Therefore, the appropriate starting point is the formal- 
ism where the cancellation of the factorial terms is taken into account auto- 
matically. This formalism is in fact none but the well-known diagrammatic 
technique for the many-body system [21], as one diagram includes the sum of 
all the processes obtained by trivial permutations. The correlations between 
different diagrams are much weaker and can be treated perturbatively. The 
subsequent sections are devoted to the statistical analysis for the many-body 
Green functions in a basis of the localized one-particle states. The transition 
resembling the Anderson's transition on the Cayley tree [17,22] will be indeed 
demonstrated, despite of the subtleties discussed above. 
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3 Choice of the model 



The purpose of this section is to introduce the simplest model describing the 
metal-insulator transition for interacting electrons as a coarse-grained version 
of the initial Hamiltonian of interacting electrons in the disorder potential. 

The Hamiltonian H of electrons placed in a disorder potential U (r) and inter- 
acting with each other via a two-body interaction potential V{r, r') — V{r' , r) 
can be written in the basis of the exact single-electron wavefunctions ^ as 



H^Ho + Vint, (43a) 
^o = E^«4ca; (43b) 

a 

^int ^\Y1 ^apjSciclc^Cs; (43c) 

Vir, r')g^s{r)Qf,,{r') d'^r d'^r'; (43d) 
gas{r)^<l>l{r)Mr), (43e) 
where cjj, and are the fermion creation and annihilation operators: 

{4; Cp} = Sa/3, {Ca, Cp} = {c),; cj} = 0, 

and {...;...} stand for the anticommutator. 

Index a labels the one-particle state and is the corresponding eigenvalue: 

.{r) = CaMr)- (43f) 



To characterize the single-particle spectrum, we introduce average density of 
states (DoS) per unit volume 

^ = ^E(^(ea)) (44) 

where (...) denote the averaging over the disorder realization and V is the 
volume of d-dimensional system. 

Both one-particle energies and interaction matrix elements are random quan- 
tities which are functionals of the random potential U{r). The usual pro- 



^ For simplicity, we consider the electrons as the spinless fermions. There is no 
reason to believe that the straightforward inclusion of the spin degrees of freedom 
affects qualitatively the final conclusions. 
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gram [15] of averaging over the disorder realizations faces difficulty when the 
relevant spatial scale becomes comparable with the one-particle localization 
length Qoc- However, we would like to discuss the physics associated with 
length scale much larger than doc- To perform this task, we adopt the reduced 
statistical model for the matrix elements and eigenfunctions which, however, 
keeps all essential physics intact. We will not discuss how the parameters of 
this low-energy effective model are connected to the properties of the system 
in the high-temperature regime [15]. 

To write down this model, we use the following properties of the localized 
non- interacting system [23]: 

(i) The typical wavefunctions have exponential envelopes, 

-ln|0„(r)| = ^^^, (45) 

sloe 

where Qoc is the localization length, and p„ characterizes the position of 
the "center of mass" of the wavefunction. ^ 

(ii) Levels repel each other for |Ca - C/?! ^ (^c^xp 

almost independent otherwise. Here we introduce the main energy scale 
of the problem: 

k ^ (46) 

which has the meaning of one-particle level spacing on the localization 
length [d is the dimensionahty of the system); 

(iii) The overlap between wavefunctions decays exponentially with the dis- 
tance, sec Eq. (45), whereas the overlap between the wavefunctions with 
the centers of mass at the distance much smaller than the localization 
length strongly depends on the corresponding one-particle energies ^a- In 
particular, for the estimate of the contribution to the interaction matrix 
elements (43d) from the distances smaller than the locahzation length 
one can approximate functions (){r) from Eq. (43e) as Gaussian variables 
with the correlation function ^ 



and are 



kl-f2| 



^ Everywhere we assume that at the distances smaller than the locahzation length 
Qoc the wave functions are well described by the semiclassical approximation, i.e., 
{mepY^'^Cioc S> 1. We will also neglect the weak energy dependence of the localiza- 
tion length as well as of the average DoS in Eq. (44). 

^ We will neglect the correlation of the wavefunctions stemming from the Cooperon 
contributions, i.e., assume the unitary ensemble 
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where cu^/? = \$,a — Cpl, diffusion lengtfi is given by — {D/u))^^'^ , and 
D is tlic classical diffusion coefficient. Dimensionless diffuson Tld{x) is 
given by 

and it decays exponentially at a: ^ 1. Equation (47) is valid provided 
\fi,2 — Pa ibI C/oc SO that the results obtained for the metallic states 
are applicable ^ . It is important to emphasize that Eq. (47) implies the 
relation between the spatial distance and the energy transfer and does 
not impose any restrictions on the values of energy themselves. 

These facts suggest the following coarse-grained version of the Hamiltonian 
(43). We discretize the space into a rf-dimensional cubic lattice with the lattice 
constant Qoc, the coordinate of each site will be labeled as p. We will call the 
unit cell of this lattice a localization cell. Each localization cell, p, contains 
large number of levels, A?" — > oo, labeled by integer Z, 

-^<«P)<^; 1<1<N. (48a) 



Two levels ii{Pi), imiPi) independent for ^ p^ and repel each other oth- 
erwise. To characterize this repulsion, let us consider the probabihty, P(n, E\ 
to find n levels in the energy interval of the width E ^ 5^. We will write this 
probability in the form 



P{n, E) 



ni 




exp 




(48b) 



The first factor in this expression characterizes the Poisson distribution of the 
independent random levels. The last exponential describes the level repulsion. 
The precise functional form of this repulsion is not important for us, we require 
only 

lim x~^V(x) = oo, 

which is a natural assumption for any repelling levels with the scale of the 
repulsion determined by (5^. 

As follows from previous discussion of the structure of the wave functions, 
the interaction matrix elements are largest for the states belonging to one 
localization cell, and decay exponentially with the distance. Thus, we will 
take into account only interaction within one localization cell. 

The interaction within one cell, however, can not cause the delocalization in 
space and, in fact, does not cause any qualitative effects at all. Therefore, we 



For d > 3 Anderson model with all the one-particle states localized the wave- 
functions at distances smaller than Q^c &re critical rather than metallic. It does not 
affect the final form of the effective model proposed in this section. 
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will take into account the single-electron hopping from one localization cell 
Pi to its nearest neighbour Being small, this hopping does not change the 
localization properties of the single electron wave function, however, taken 
together with the electron-electron interaction would eventually lead to the 
transiton. 

The resulting Hamiltonian, thus, takes the form 



H = Ho + Vint; (48c 



(48d) 



^0 = E 4ip) 

p,l 

Vint ^ I E ^/r(p)4(p)4(P)^..(P)2.i(P)> (48e) 

hl2jij2;p 

where {cl{pi); Cjip^)} = SijSp^p^, {q(Pi); 9(^2)} = {cl(pi); c](P2)} = 0- 

Here a arc the vectors connecting the cell p to its nearest neighbors. Dimen- 
sionless hopping parameter, /, such that 



is introduced to control further perturbative expansion. We will chose / > 0: 
this choice does not affect any conclusions, as / will connect terms with the 
random signs. 

The antisymmetrized coefficients Vi'j^^{p) = Vl^i^^{p) = — ^^^^^^^(p) are random 
numbers. Because the physical processes discussed in this paper are associated 
with the counting of the resonant denominators, the particular choice of the 
statisitical distribution of Vi^il^{p) is not really important. For the calcula- 
tional convenience we choose the binary distribution 



where A ~ / -C 1 is a dimensionless parameter allowing to control the pertur- 
bative expansion, ujij — $,1 — and we omitted argument p on both sides of 
the equation. We chose A > without loss of the generahty, as all the effects 
which will be considered are not sensitive to the sign of A. Function T{x) is 
introduced to describe the interaction decaying rapidly with the distance be- 
tween the levels in the energy space, see Eq. (47). As the smallest linear scale 
in the reduced model is Qoc the maximal energy transfer which is permissible 
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to consider in the model is of the order of 5^. Thus, we choose 




(48h) 



Actual value of the parameter M is not well defined and we will consider it 
as the initial data for the effective Hamiltonian. The significance of the upper 
bound to M for the consistency of the theory will be clear later, see discussion 
after Eq. (80) as well as Sec. 7.1. 

The signs for different wave functions are not correlated, so 



Equations (48) constitute a complete formulation of the reduced model for the 
interacting electrons in system with localized one-particle states. This model 
will be analyzed in the subsequent sections to show the stability of both high- 
temperature phase (metal) and low-temperature phase (insulator). 



4 Formalism. 

The purpose of this section is to describe the machinery which enables us to 
put the previous qualitative arguments into the context of the usual many- 
body theory of non-equlibrium systems. We will start with the outline of the 
Keldysh formalism for the exact (non-averaged) Green functions correspond- 
ing to Hamiltonian (48) in Sec. 4.1 and formulate which quantity describes 
the metal-insulator transition in Sec. 4.2. Next, in Sec. 4.3 we will describe 
our main working approximation which corresponds to the summation of all 
the rainbow diagrams (SCBA). The justification of the validity of this approx- 
imation for the description of the transition is postponed until Sec. 7. 

4-1 Time evolution equations and basic definitions. 

We intend to describe both metaUic and insulating regimes. In the latter 
regime relaxation dynamics is absent, there is no mechanism to establish the 
thermal equilibrium, and the temperature itself is not defined. Therefore, the 

^ It can be shown that taking into account the algebraic decay of T(x), a; S> 1 
would require the consideration of the spatial correlation of the wave functions, so 
the approximation for the interaction matrix elements to be independent of each 
other and of the hopping would be false even on a qualitative level. 




1; {(^i{p)4 ip')) = Sp,Ai'^n'- 



(48i) 
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only appropriate formal framework is the non-equilibrium (Keldysh) formal- 
ism [24]. We define the corresponding Green functions as 



gp{h,t2;p) = -ie{ti-t2)(^{ci{h,py,cl{t2,p)}J^, 

g^ih, t2; p) = ie{t2 - p); c\{t2, p)} )), (49) 



g^{hM\p)^-i{{ ci{h,py,cl{t2,p) 



where 6{t) is the Heaviside step function, and the fermionic operators are writ- 
ten in the Heisenberg representation. Quantum mechanical averaging {{■■■)) 
is performed over an arbitrary density matrix to be found from the solution 
of the kinetic equation. To avoid misunderstanding, we emphasize that no 
averaging over the disorder realization is assumed in Eq. (49). 

We parametrize the Keldysh Green function as 

gl'ip) = gp{p)oni{p) - ni{p)og^{p) (50) 

where we omitted the time arguments for brevity and introduced the short- 
hand notation 

CoV = I dhCih,ts)V{ts,t2). (51) 
for arbitrary functions C, T>. 

In the thermodynamic equilibrium 

ni{p;e,t)^l-2Me), (52) 

where /^(e) is the Fermi distribution function with arbitray temperature and 
chemical potential, and the time Wigner transform is defined as usual: 

In the absence of interaction ni{p; e = ^/(p, t)) characterizes the occupation of 
the level ni{p; ^i{p),t) = 1(— 1) for an empty (filled) level (p, /). 

In what follows we will use standard diagrammatic technique for the pertur- 
bative expansion for the model (48). The basic elements of this technique are 
defined in Fig. 3. 

The Green functions (49) corresponding to the Hamiltonian (48) satisfy the 
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i \Gi{p)] 



Gf(p) Gf (p)' 
Gf(p), 



{l',p')^il,p) 



Gf (p) = [Gf (p)l = 



■6(p) + i0+ 



Gf{p) = -2Trini(p,e)5[e-^i{p)] 



b) 



Ml 



c) 

(^l>Pl) I {l2,P2) 
' Ml 



d) 




e) 





il:i(pj = Fig. 4 +Fig.5 + ... 



Fig. 3. (a-d) Basic elements of the diagrammatic technique. The Keldysh space is 
labeled by iJ,i, and the Pauli matrices in Keldysh space are defined in Eq. (54). 
Line crossing the Green function excludes the orbital (I, p) from the summation, (e) 

Representation of Eqs. (54) - (55). 

equations (we omitted time arguments for brevity) 

[idt, - 6(P)] Qi{p) = r^Kti - + ±i{p)ogi{p)- 
[-idt, - 6(p)] Qi{p) = To5(ii - t2) + gi{p)o±i{p); 



Qf{p) 



-I K 



^lip) (p) 

5]f(p) 



K 



1 




1 








1 




1 




K 





K 



(54) 
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Fig. 4. First (a), second (b) and third (c) orders contributions to the self-energy 
S;(p). The diagrams which will be taken into account in the self-consistent Born 
approximation are highlighted by the dashed frames. Cross-sections of the diagram 
which produce the imaginary part of the self-energy are denoted by the dotted lines. 

The self-energy ^^(p) is given by 

±i{p)=fj + ±i{p), (55) 

where fj originates from the coupling of the sysstem to an external bath with 
regular continuous spectrum (it will be discussed in more details later). This 
coupling has to be kept small but finite and can be put to zero only in the 
end of the calculation. Self-enery S;(p) originates from the electron-electron 
interaction and hopping and represents the sum of all diagrams which cannot 
be separated by cutting one-electron line (/, p), see Figs. 4, 5. 

Substitution of Eq. (50) into the Keldysh component of Eq. (54) yields 

{dt, + dt,) ni{p) = -iSf (p) o ni{p) + ini{p) o ^f{p) + i^f{p). (56) 

On the next step we restrict ourselves to the consideration of very slow dynam- 
ics. In this case one can perform the time Wigner transform (53) in Eq. (56) 
and obtain the quantum Boltzmann equation 

dtni{e- p, t) + [S^iie; p; {n}) * ni{p, e, t)] = St, (e; p; {n}, {dtu}) + St|'«*^ (57) 
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cl c2 c3 cl c2 c3 



Fig. 5. Fourth-order contributions to the self-energy S/(p). The diagrams which 
will be taken into account in the self-consistent Born approximation are highlighted 
by the dashed frames. Cross-sections of the diagram which produce the imaginary 
part of the self-energy are denoted by the dotted lines and denoted by cl - c4. The 
crosssection c4 implies cutting three lines at the same time. 

where the collision integrals are defined as 

St,(e; p; {n}, {dtfi}) = -2n(e; p; {n}, {dtn})ni{p, e, t) + iSf (p; e; {n}, {dtfi}); 
^^bath ^ _2 ini r)^ni {p, e,t) + irj^ . 

(58) 
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Hereinafter, the time Poisson brackets are defined as 

[C*V] = dtC dj) - d,C dtV. (59) 

for arbitrary functions C{e,t), V{e,t). The entries in Eq. (57) are defined by 

T.^ = 5i- iV- = 5^ + iV- (60) 

where we suppressed all arguments which are the same as in Eq. (58). Using 
Eqs. (60) and the analytic properties of the retarded and advanced Green 
functions one obtains 



g^{p-e,t) = [gt{p-e,t) 



de'Aiip- e',t) 



e — e 



i0^ 



(61) 



TT ^Ti{e\p]{n},{dtn}) 



[e - Up) - p; {n})Y + [F^e; p; {n}, {dtu})] 



2 ' 



and = Ti-\-lm.r]f. Notation T{{nY) means that T is b, functional depending 
upon all the functions n;(e, p) but local in time. The latter functions enter the 
expressions through the quasistationary version of Eq. (50): 

gl'ip;e,t) = -2mni{p;e,t)Ai{p;e',t)-z J de' [Ai{p;e' ,t) * ni{p;e,t)] P-^, 

' (62) 

where P denotes the principal value. 

What remains now, is to specify the thermal bath. As we already mentioned, 
the particular form of this choice is not important. We will require only that 
it preserves the number of particles 



J deAiip; e,t) St''''\e,t) = 0, 



it is local, and the collision integral St^"*'* is nullified by the equlibrium distri- 
bution function (52). It is easy to check that the choice 



2Im7^'^ = j dujujh{uj) Ai{p;€ — u;,t) 
dujujb{uj) Ai{p; e — uj,t) 



UJ 



coth — + ni{p;e-u;,t) 



-IT] 



UJ 



coth — ni{p;e-uj,t) + 1 



(63) 



where b{u;) — b{—u;) > 0, saisfies both those requirements^ . 



^ The subsequent formulas for non-zero coupling fj actually may be used to describe 
the effect of the short-rangle electron-electron interaction on the phonon-assisted 
nearest neighbors hopping. 
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4-2 Formulation of the problem. 



Having introduced the definitions, we are ready to reformulate the criterion 
distinguishing insulating and mctalUc states. The left-hand side of the Boltz- 
mann equation (57) describes the evolution of the occupation of the levels in a 
self-consistent field created by all the other electrons. This evolution is deter- 
ministic and time reversible. It is the right-hand side of the kinetic equation 
(collision integral) that makes the time evolution probabilistic and specifies 
the direction of the time arrow. Thus, the energy dependence of the decay 
rate 2ri{e) determines whether or not the irreversible evolution occurs in the 
system. If coupling with the environment, b{uj) from Eq. (63), is finite, Ti{e) is 
positive for any energies. However, if this coupling tends to zero (but not faster 
than the exponential function of the volume of the system V, 6 > exp [— V/Vq]), 
two situations are possible, see Fig. 6a,b. 

(i) Despite 6(a;) tending to zero the number of the intermediate states via 
which the excitation can decay goes to infinity. This results in r;(e) 
being a smooth function of energy even at 6 (a;) — > 0. This situation 
corresponds to the applicability of the Fermi golden rule, the thermal 
equilibrium within the system established at times independent of the 
external bath, so it is natural to classify this regime as metallic. 

(ii) The number of the intermediate states via which the excitation can decay 
remains finite and independent on b{u!) as 6(a;) — > 0. This results in 
r;(e) to be a sequence of resonant peaks positioned at the energies of 
the exact excitations of the many-body system. In this case n(e) will 
remain extremely singular function whose relaxation rate is determined 
by 6(a;). At 6(a;) the thermal equilibrium can never be reached and 
this regime is insulating. 

Since these two behaviors arc qualitatively different, there could be no smooth 
crossover between them, and only phase transition is possible. Therefore, to 
show the existence of the transition, it is sufficient to formulate the conditions 
at which either metallic (i) or insulating (ii) regimes are stable. It will be done 
in Sees. 5 and 6. The investigation of the behavior of the kinetic coefficients 
near the phase transition point itself will be the subject of a separate paper. 

Next question one has to ask is how to distinguish between metallic and insu- 
lating phases within a statistical framework. It is clear that the positions of the 
peaks in r(e) for the insulating regime deviate randomly with the variation of 
random energies from Eq. (48). Therefore, the averaged value of the decay 
rate (r(e)) is qualitatively similar in both phases and can not be used for the 
distinction. The magnitude of the fiuctuations ([5r(e)]^) is, however, qualita- 

Subsequent discussion is a straightforward generalization of the Anderson's argu- 
ment [1] to a many-body system. 
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Fig. 6. Schematic energy dependence of the quasiparticle decay rate r(e) for (a) 
metallic and (b) insulating phases. The corresponding distribution functions are 
sketched on panels (c), (d). The characteristic functions are plotted on panels 
(e), (f). 

tively different in both cases. In the metallic case the averaging is performed 
with respect to the smooth positive functions, whereas in the insulating regime 
the fluctuations are determined by the squares of the separated delta-peaks 
and thus diverge as the width of the delta peaks goes to zero. Thus, we have 
the criterion 



([^r(6)]^) 



lim lim 



finite; metal 
oo; insulator. 



(64) 



Another way to address the same problem is to investigate the distribution 
function P{r), see Fig. 6(c-d). One finds by simple inspection of Fig. 6(a-b) 



lim lim P(T > 0) 

&(a;)-»0 V-»oo 



> 0; metal 
0; insulator. 



(65) 
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For us, it will be more technically convenient to perform the Laplace transform 
and calculate the characteristic function 



W{s) = {cxp [-sT{e)]), 



(66) 



where the precise definition of the averaging procedure is given in Sec. 6.1. 
The criterion (65) [see also Fig. 6(e-f)] of the insulating phase translates into 



for any fixed s > 0. 

Closing this subsection, we emphasize that the transition occurs as a function 
of temperature (which corresponds to the extensive energy of the many-body 

state) and not as a function of e which charaterizes the energy of the one- 
particle excitation on top of this many-body state. The latter energy is not 
an extensive quantity and can not be a characteristic of any phase transition, 
in contrast with conclusions of Ref. [16] for a finite-size system. 

4-3 Self- Consistent Born Approximations (SCBA and ImSCBA). 

In this subsection we introduce our main approximation for the sTimmation 
of the infinite series of perturbative expansion. We will discuss motivation for 
this approximation here in quite loose terms and justify it further in Sec. 7. 

Contributions to the self-energy shown in Figs. 4,5 bear different physical 
significance. For instance, the Hartrcc-Fock diagrams (b), (d2) of Fig. 4, and 
diagrams (3) of Fig. 5 characterize the self-consistent one-particle spectrum. 
As we explained before, the structure of this spectrum is not relevant for 
the transition. On the contrary, diagrams (cl) and (c3) of Fig. 4 can lead 
to irreversible processes, and the appearance of the imaginary part of those 
diagrams signals the metal-insulator transition. Thus, those diagrams will be 
taken into account. 

The third order diagrams. Fig. 4d, describe the effect of the change of the self- 
consistent potential on the tunneling process (dl); effect of the tTinncIing on 
the self-consitent potential (d3); and the effect of the self-consistent potential 
on the spectrum of the decay channel. These diagrams will be neglected for 
the reasons explained above. 

Fourth order diagrams Fig. 5 describe further potentially irreversible process: 
tunneling out of the localization cell (1); three-particle production with the 
consequent tunneling (3); five-particle production (5). Those contributions 
must be taken into account. On the other hand, diagrams (2) and (4) once 



lim lim Wis) — 1 



(67) 
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(c) 



{h,Pl) {l2,P2) {h,Pz) 

{h,Pi) / {h,p3); {l2,P2) / {l,p) 




Fig. 7. a) Self-consistent Born approximation. The thick fermionic hnes are defined 
in Fig. 4a. b) Spurious contributions of the fourth order generated by iteration of 
SCBA: b.l contradicts the definition of the self-energy; b.2 contains the intermediate 
particles in the same states; such contributions are cancelled by diagrams (5.1a) and 
(5.2a)-(5.2c) of Fig. 4. c) Fourth order tunneling and interaction contributions not 
taken into account by SCBA. These contributions have a random sign. 
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again represent the self-consistent potential affecting the lower order processes 
(cl), (c3) of Fig 4. Diagrams (6.1), (6.2) of Fig. 5 describe the interaction of 
the quasiparticles created by the lower-order process (c3) of Fig. 4. Such effect 
has the same physical significance and the relative value as the effect of the 
self-consistent potential and will be neglected. The effect of those diagrams 
will be estimated in Sec. 7. Finally, diagrams (5.1a), and (5.2a)-(5.2c) are 
exchange counterparts of the main diagrams (5.1)- (5.2) [we will elaborate it 
further in Sec. 7]. Their role is to cancel out contributions in (5.1), (5.2) which 
are forbidden by Pauli principle (more than one particle in one intermediate 
state) and also to give the random sign interference corrections. The latter 
corrections will be neglected here and estimated further in Sec. 7. 

The above discussion of the lowest orders of the perturbation theory suggests 
the following prescription for the summation of the leading series: (i) we take 
into account only the even orders of the perturbation theory; (ii) we require 
that each contribution maximize the sum of the tunneling events and the 
extra quasiparticle production; (iii) we neglect contributions with random sign. 
By inspection, one can see that such series with the correct combinatorial 
coefficients is generated by the self-consistent Born approximation (SCBA) 
shown on Fig. 7a. 

Iterations of the SCBA equations also produce spurious contributions shown 
on Fig. 7b. For the SCBA scheme to be vahd, we will make sure that those 
spurious contributions are always smaller than those responsible for the final 
result. 

We calculate the self-energy (Fig. 7a) according to the rules of Fig. 3. We 
neglect the shifts 5^i{p) in accord with our previous discussion, i.e., we make 
the additional approximation to the usual SCBA scheme 



Er"")(p;e)-iImEr-^)(p;e) (68) 



We will refer to approximation (68) as Im-Self-Consistent-Born-Approxima- 
tion (ImSCBA). In the two subsequent sections, we will work only within this 
approximation. 

Using spectral representation (61) and (62) and performing the integration 
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over the intermediate energies, we obtain 



(in). 



^(bath) I 



h,a 



rf")(e) = 7rA2<5^2 ^ y^j^ J ^e, de, de, A,ie,) A^e^) M^s) 

X S{e -61-62 + 63) F^i^.i^{6i, 62\ 63); 



rr'^(e) = ^/ dujuh{u)Ai{6-u) 



coth — +ni(e-u;) 



(69a) 



[e-6r + [rz(e)]^ 



T 



62 ~ 



T 



T 



61 - 



T 



■^^2 - 63 



1 2 



^iti2;/3(^i'^2;e3) = ^^1 + ni,{6i)ni^{62) - rii^iea) K(ei) + n;2(e2)] |; 

where 77 is defined in Eq. (63), and we utihzed the notation used in Eqs. (48) 
for the nearest neighbours. Everywhere, the coordinate p and time t are as- 
sumed to be same in all terms in the equations unless it is specified explicitly 
otherwise. 

Equations (69a) form a closed set for finding the decay rate for fixed occupa- 
tion numbers ni{e,p). In the delocalized regime the time evolution of those 
occupation numbers is governed by the ImSCBA version of the kinetic eqution 
(57). Calculating the Keldysh component of Fig. 7a and using Eqs. (58) and 
(63), we find 

dtni{6)^Stf^ + Stt^+Stf''"'^; 

Stf^ = 27rP5l Yl (e, P + a) (e, p + a)-ni (e, p)] ; 



ll ,a 



Stf"^ = 2t:X'5. 



2x2 



H ^h'L j ^^1^^2d6zAi^{6i)Ai^{62)Ai^{63)5{6-6i-62 + 63) 

^it/2;i3(^i'^2;e3) = ^1 - nicies) [1 + ni,{6i)ni^{62)] + K(ei) + ^^^(es)] |; 

^^{bath) ^ I d^^^(^) Ai{e-Lo) 

x jcoth ^ [ni{6 -uj)- ni{6)] + 1 - ni{6 - uj)ni{6)^ . 

(69b) 

Equation (69b) is the usual quantum Boltzmann equation written in terms of 
the exact (not averaged) single electron levels Ci{p). As any Boltzmann equa- 
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tion, it must respect the fundamental symmetries of the system: conservation 
of the number of particles for any collision; conservation of electron energy for 
the processes involving the electrons only; and the conservation of the number 
of particles for a given energy for elastic coUisions. It is straightforward to 
check that the collision integral indeed possesses the desired properties: 

i,P 

J2 J deMe,p)Stt\e,p) = 0; 
dee^Ke,P)St^He,p) = 0; 

I •' 

j:[deMe,p)Stf'^"'\e,p)^Q. 
I 



The properties (70) of the collision integrals enable one to write the continuity 
equations for the particle and energy densities H and and introduce the 
corresponding currents J, Jg: 



dtM{p) + dwJ{p) = Q; 

dtS{p) +d\Y J s{p) = -Qbathip); 



(71a) 



where Qbath is the thermal flow to the thermal bath; it will not be important for 
the further cosideration. The lattice version of the divergence of the currents 
is defined as 



1 



divJ(p) = — ^ 

<5«OC ji. = l 



a 



a 



(fc)^ 



and it becomes the usual divergence in the continuum limit. Index k labels 
the direction in the cartesian coordianat system in d dimensions. Vector a^''^ 
is the lattice vector along the kth direction. 



The densities are defined on the sites p 



s 



de 



EA(6,p) 



1 -n/(e, p) 



(71b) 



whereas the currents are defined on the links p*^*^' = P + 



jk 



de 



X 



^loc 



J2 Me, p) Ai,(e, p + a^'^) 

l,h 

ni,[e, p + a^''^) - ni{e, p) 



(71c) 
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Finally, St}^ is nulled by any function n/(e; p) = n(e), electron-electron inelas- 
tic collision integral St^*"^ is nulled by any Fermi finction ni{e; p) = tanh 

and the St;^"*^"* vanishes for n;(e; p) = tanh^^^. 

The closed system of Eqs. (69) is a drastic (though parametrically justifiable) 
simplification in comparison with the original problem. However, it is still 
a non-linear system which depends on an infinite number of random energies 
ii{p). Substantial progress can be achieved within the statistical analysis. This 
analysis is a subject of two following sections. 



5 Stability and properties of the metallic pheise. 

5. 1 Condition for stability 

The hallmark of the developed metallic phase is self-averaging of the kinetic 
coefficients, see Sec. 4.2. To establish a sufficient condition for the existence 
and stability of this phase, we assume that the inelastic decay rate is indeed 
self- averaging and then justify this assumption by explicit calculation of its 
mesoscopic fluctuations, see Fig. 8. As a result, we will see that this condi- 
tion reduces to a certain integral inequality for the distribution function, see 
Eq. (77). We will also see that even when this condition is satisfied, the elas- 
tic rate still may not be self-averaging. This, however, will not violate the 
criterion (64), as the fluctuations of the elastic rate remain finite. 

First, we take into account only the inelastic rate in the Ai of Eq. (69a). We 
assume and justify a posteriori that the main contribution originates from 



Statistical averaging over the distribution of the levels can be performed 
independently. 



For the needed products of the spectral densities, see Eq. (69a), we find 



I - k\, \k - lj\ > 1, i,j = 1,2,3. 



(72) 





(73) 



1 2r(^")(ei) 



^'^c(ei-e2)2 + [2rM(ei)] 



2 ' 
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Fig. 8. a) Simplification of tlic self-consistent Born approximation at T > Tj„. b) 
The self-averaging of the hopping term at T > T^i. Cross with the dashed line(s) 
stands for the averaging of the corresponding Green function (product of the Green 
functions) over the realization of random energies ^i{p), see Eq. (48a). 

where we assumed that r(e) is a smooth function on the scale of r(0). 

Averaging Eq. (69a) with the help of Eq. (73), sec also Fig. 8a, and assuming 
that the distribution functions ni{e, p) do not depend explicitly on the orbital 
index /, we find: 



(rf")(6)) ~ / deida;F^(6i+c.,6i;6-a;) 



X 



ei + u - e 



(74) 



As we will see shortly, the metallic regime is realized when the characteristic 
energy scale of the distribution function is much larger than M5(^ [Eq. (48h)] - 
typical transferred energy in Eq. (74). Under such conditions, the second term 
in the second line of this equation is more restrictive on the phase volume and 
that is why it can be neglected: 



(r^He)) 



deida;T^(^^ 1 F= 



(ei +a;,ei;e - a;). 



(75) 



Calculating mesoscopic fluctuations of the inelastic rate shown in Fig. 8b with 
the help of Eq. (73), and keeping the terms with the largest phase volume we 



^■•^ All the formulas of Sec. 5 are written under the condition (r^^'"'^(e)) <C M5(^. In 
this case one can neglect the dependence on the index I of r^(e), ni^e), and the 
collision integral. Such dependence is present only for |e — Cz| ^ ^-^<^C' while all 
physical properties are determined by the region particular, the 

distribution function enters the observables only as J2i ^li^) ^i(^)- 
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find 



(in), 



I 



dei do; 



r{m)(e + 00) + rM(ei) + W{ei - u) ' 

(76) 



Our initial assumption that the inelastic rate is self-avergaing is justified pro- 
vided that 



< 



^(m) f 



(77) 



According to Eqs. (75) and Eq. (76), both sides of this inequality are deter- 
mined by the distribution function only, so that Eq. (77) is a sufficient (but 
not necessary) condition for the arbitrary non-equilibrium state to be metallic. 

For the thermal distribution n(e) = tanh^, the explicit expressions can be 
obtained. One finds 

^rM^ = nX^MT, (78) 
where M is the coefficient defined in Eq. (48h), and 

The condition (77), thus, reduces to the lower bound for the temperature 

To complete our discussion of the properties of inelastic rate, we justify our 
assumptions. To check the condition (72). we anlalyze the structure of the 
energy integrals in Eqs. (74)-(76) and find |/ — I/2 — ^sl — ^ S> 1, |^ — /3I — 
T jbc^ » M, which is consistent with Eq. (72). Deriving Eq. (78) we assumed 
T » Mbc^. It is consitent with Eq. (80) provided that the condition (48h) is 
fuUfiUed. 

Let us turn now to the properties of the elastic decay rate, r*^^'\ of Eq. (69a) 

r;^')(e,p) = ^<^?^'E^(^'P + «); 

a 

, 1^ r(^")(e) (81) 

^ r - ejp)f + [r(^-He)f 

Let us note that only the inelastic width enters the right-hand sides of these 
equations. Further iterations of the elastic processes in expression for A/(p+a), 
see Eq. (69a), generate either terms small as /^(r*^*")(e)/5^) -C 1 (originating 
from off-resonant levels), or the divergent term corresponding to the elastic 
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^(P2;e) 



Sn{p-^;e) A{p2;e) Sn{p^;e) 



Fig. 9. Sketches of the shapes of the non-equihbrium influx function 6n{pi,e) and 

the spectral density of the other neighbours A{p2;€) = X^;^i(P2)^) of the other 
neighbour; for a) low-temperature metal Tin < T < T^i; h) T > T^i. See text for the 
further explanation. 



return on the same level. The latter contribution, however, corresponds to the 
spurious diagram Fig. 7b. 1 and must be discarded. 

Average and fluctuations of the elastic rate (81) is calculated with the help of 
the Eq. (73) and we flnd 



2\ _ dTrI^6l 



r(m)(e)' 



(82) 



where 2d is the number of the nearest neighbors. 



Using Eq. (78), we find 



el 



167r2dMA2 A 



(83) 



where the numerical factor is chosen for the convenience in the further formu- 
las. At T » Tgi the level discreteness plays no role, so that both elastic and 
inelastic decay rates arc self- averaging, see also Fig. 9. At Tin <^Tei only 
p(«n) jg self-avcraging, while the fluctuations of F*^"^'^ are large compared to the 
average. However, they are flnite, so that the system is in the metallic state 
according to criterion (64). 
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5.2 Kinetic equation and transport coefficients. 



Having studied the statistical distribution of the decay rates, we are ready to 
apply the same ideas to the kinetic equation: 



n{e, p + a)- n{e, p) 



(84) 



9tn(e, p) = 27r/\^ 5^ A(e, p + a) 

a 

+ ^/du;d6iT^(^) A(6 + a;,p) 

X I [n(e + a;, p) - n(e, p)] [1 - n(ei, p) n(ei - a;, p)] 
+ [1 - n(e, p) n(e + p)] [n(ei - a;, p) - n(ei, p)] | 



Equation (84) enables us to find the kinetic coefficients in the system. We look 
for the distribution function in the form 



n(e, p, t) = tanh ^ + $(e, p, i) + (/?(e, p, i), 



(85) 



where the function $(e, p, t) describes the shape of the distribution function on 
the energy scale e > M5^, whereas the function (/^(e, p) encodes the structure 
on the scale e ~ Namely, we impose the condition 



e+A 



de 

{^{^,p))e= j ^<^(ei,p) =0, 



(86) 



e-A 



where A is an energy interval 5^ <^ A < M5^. We substitute Eq. (85) into 
Eqs. (84) and linearize with respect to $ and </?. For the smooth part of the 
distribution function we find 



a,$(e,p) = 27r/\^ 

a 

+ 27r7\^$:(^(e, p + a) 
+ St$$(e,p) + (St^<^(e,p) 



$(e,p + a) -$(e, p) 

^(e,p + a) - (/?(e, p) 



(87a) 



Deriving Eq. (87a) we used the fact that {A{e))^ = 1/5^ is a self-averaging 
quantity. 

The smooth part, $, contains, in particular, density and energy density which 
propagate diffusively through the system. In contrast, the oscillatory contribu- 
tion decays due to the inelastic processes. Thus, function (/? can be considered 
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in the stationary limit 



0=27r/\2^P(6,p + a) 



A{e,p + a) 



$(e,p + a) -$(e, p 
(p{e,p + a) -f{e, p) 



(87b) 



+ 



St^(^(e, p)\ 



where we introduced the notation for the oscillatory part of the expression 

II ■ ■ ■ I lose (■ ■ ■ )f ■ 

The smooth part of the linearized collision integral St<^(^ is given by 



St^(^(e))^ = 27rA2r j dwT^ (^^ j {A{e + uj) [ip{e + uj) - (^(e)])^ 
«.2r(^")5c(>l(6)(^(6)),. 



(87c) 



Here we used the fact that {A{e + u})ip{t)) ^ = {A{e + Lj))^{ip{e))^ = 0, for 
I a; I > and the contribution to the integral is determined by ~ M5^. We 
also used the expression for the inelastic rate (78). By the same token, we find 
the oscillatory part as 



j du;T^(^^^ [(^(e + u;)-(^(e)]«-2r(*'^V(e)- (87d) 



We remind the reader that the relationship T ^ Tj„ ^ M5c_ is widely used, see 
also Eq. (80). Finally, the shape of the smooth distriution function is stabilized 
by the linearized collision integral 



St$$(e) 



27rA2 



UJ 

+ 2 



dc^T^(^j{|coth^[$(6 + a;)-$(6)] 
$(e + u) tanh ^ + $(e) tanh 



If e — UJ e 

H — tanh tanh — 

4 V 2T 2T 



X j dei 



$(ei + uj) tanh ^ + $(ei) tanh 



(87e) 



This hnear operator has two zero modes 



$^(e) = -d, tanh ^; $r(e) = ^ytanh ^, 



(87f) 
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which reflect the conservation of number of particles and energy by the inelas- 
tic processes. 



Because of the condition T ^ M(5^, Eq. (87e) can be rewritten in the energy 
diffusion approximation 



St*$(e) = L'(^)<j9,^$(e) + | 



$(e) tanh — 



- — tanh [ dei$(ei) tanh > ; 
4 2TJ ^ ^ 2rj ' 



(87g) 



D^'^ = ^ . 

12 



The simple form of the collision integral (87d) enables us to solve Eq. (87b). 
We adopt the following relation 



(88) 



which is automatically fuUfiUed^^ ^j. ^ rp(in)^ ^gg^^ 7 ~ A. 

Then r(^") > PSlA{e) and Eq. (87b) yields 



X 



$(e,p + a) -$(e, p) 



l + O 



(89) 



The regime F^*") <C 15^, which may occur for somewhat artificial for our model 
choice 7 S> A, corresponds to coherent oscillations of the population between two 
resonant levels. In this regime one should modify Eqs. (89) and (90) as 

a 

2TTp5lA{e, p + a) 



2rM + 2TTPS'^[A{e, p + a)+ A{e, p)] 
dt^ie, p) - St$^>(e, p) = Y, [^(f ' P + «) - ^(e. P) 

a 

I 27r72<52^(e,p)2r(*")A(e,p + a) \ 
^ \ 2rM + 2i,P5l [A{e, p + a)+A{e, p)] / 



One should also use the exact spectral densities for the resonant pairs, as obtained 
from the solution of a 2 x 2 problem. The results for the low- and high-temperature 
regimes remain unchanged, but an intermediate region appears. This regime is anal- 
ogous to that studied in Ref. [25] for the conductivity due to the scattering of 
localized electrons on phonons. 
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Substituting Eq. (89) into Eqs. (87c) and (87a), and keeping once again only 
the terms leading in (/(5(^/r(*"^)^, we find the equation for the smooth part of 
the distribution function 

dMp) = St$$(p) + 27rI^Sl (^(e, P)^(e, P + a)), [Hp + a) - $(p) 

(90) 

where we omitted the energy argument of the distribution function. 

Equation (90) deserves some discussion. The first term on the right-hand side 
determines dynamics of the population at different energies e on the site p. The 
second term describes transport between neighboring sites. The inspection of 
this term shows that for F*^*'^^ <C 5^ this term is contributed mostly by rare 
overlaps of the peaks in the spectral densities on the neighboring sites, see 
Fig. 9. We will call such rare events "pin-holes". The typical energy separation 
between the pin- holes is 5^/F(^"). Thus, for a given energy this term may not 
be self-averaging, and the calculation of the observable transport coefficients 
requires further statistical analysis. We note that the situation when transport 
is dominated by rare configurations is quite common in systems with strongly 
fluctuating local transmission, see Ref. 



To perform this analysis, we assume that the energy relaxation rate within 
each site is larger than the rate of the tunnelling into the neighbors (this 
assumption is justified for our model in Appendix A). In this case the solution 
is restricted to the zero modes of the inelastic collision integral (87f) 

$(e, p; t) = 5i^{p, t) $^(e) + 5T{p, t) ^^(e). (91) 

Substituting Eq. (91) into Eq. (90), performing the corresponding energy in- 
tegration to utilize the properties / deSt$$(e, p) — and / deeSt$$(e, p) = 
0, and neglecting thermopower and Peltier coefficients (which have random 
signs), we find 

e^dfi{p) , a 



Sc dt 
dT{p) 



Eg(p + ?) [Mp + «)-Mp)]; 

" (92a) 
^ = E^(P + f ) [Tip + a) - np)] ; 

where Cy = Tr^/yT/S is the specific heat per localization cell. The electrical 
conductance G and the thermal conductance /C are defined for each link p+a/2 
connecting two sites p and p-\- a. They are given by 



a 



(92b) 
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where Ba,K are dimensionless random quantities determined by the overlaps of 
the densities of states: 

ct\ „n /■ de 



^^S^g (W\ 2^,rM/vr (92c) 

jjf2T^'''\ 2T ) [6(p)-6'(p')]' + [2r(-f' 

where p' = p + a. Wc used the explicit expression Eq. (81) for A{e^ p) and the 
condition V^'''^^ -C T. Dimensionless functions Pa,K{x) such that / Pa,K{x) dx — 
1, are given by 

2 cosh X TT^cosh x 



Equations (92) are nothing but the description of a network of random con- 
ductors (thermal conductors). It is easy to check that the average of these 

quantities over the realizations of $^1 gives the temperature-independent result 
(Bct^k) = 1- However, this result is meaningless as the observable conductiv- 
ities are determined by typical rather than rare events. In fact, observable 
conductivity a and observable thermal conductivity k are given by 



h 



X < 



exp[(ln,B^) 



d = 1; 



d = 2; 



exp 



(ln"^-B<,)V"<il ■d>2- 



d= 1; 



exp[(lni3,)]; d = 2; 



exp 



{In'''' B^y/"A ]d>2; 



(93) 



where d is the dimensionality of the system and are the numerical coeffi- 
cients of the order of unity which are not known analytically for d > 2. The 
formula for 0? = 1 is the trivial result for the random resistors connected in 
series whereas the result for 0? = 2 follows from the duality arguments [27] . 



The averages entering Eq. (93) can be immediately calculated if the charac- 
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teristic functions Pa,K{s) — (e ^^'^■'^N of the distributions are known: 



B. 



dsP„,,{s); (lni3,,,) = / — e"^ -P.,.(s) . (94) 



The characteristic functions Pcr,K{s) can be found from the definition (92c) by 
a straightforward calculation given in Appendix B: 



PaA^) = exp i - 



dx 



rSo 



, 52(y) = ye-^[To(y)+Xi(y)], 

(95) 



where Tq and Ti are the modified Bessel functions, and (5cr,K are given by 
Eq. (92d). The result is controlled by a dimensionless parameter 



r(T) = 



(96) 



The meaning of this parameter is the typical number of resonances in the 
energy strip |e| < T". 

The result of numerical integration of Eqs. (94), (95) is plotted on Fig. 10. In 
two limiting cases r ^ 1 and r <C 1 we have 



r > 1 : { 



r < 1 : 



/V(s)f^e-^[l + (l/6)sVr], 

P^{s) ^ e-'[l + (7/10 - 6/7r2)s7r], 

PJs) ^ e-^. 



(97) 



where G ^ 0.916 • • • is the Catalan's constant. Substituting Eq. (97) into 
Eq. (94) we find that for r » 1, Pa,K{B) are strongly peaked around 1, so that 



13 



The average of the logarithm is calculated using the identity 



/ 



ds 



In 
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the temperature indepedent Drude resistivity is restored: 





24 \ 








rp2 



Too- ^ , «:oo(Jj- 3^ 



(98) 



In other words, the effect of localization of one-particle wave functions on 
transport is completely removed by the inelastic processes even in the "non- 
ergodic" regime of F'-*"^ <^ 5^, where the peaks in the one-particle density of 
states are still well-resolved. 

In the opposite case, r <S 1 we have for d = 1, 2: 

^2 



TT / TM nr^^r. 48G' / T' 



a{T < ^JS^T^l) = (Too- ( jjr j, i^{T <^ \JkTei) = i^oc{T) {ST T 

\ C / \ C / 

(99) 

(for larger dimensionalities the temperature dependence is the same but the 
numerical coefficients could not be found analytically). 



These results have the following meaning. At T ^ Jd/T^i the quantities are 



self- averaging and the result is temperature independent. At T ^ ^5cTei the 
electron on a site cannot explore enough pin-holes to find the rare resonant one 
which would determine the average. As a result, it chooses the best available 
pin-hole {i.e., the one with the smallest separation between the levels). The 
denominator in the Eq. (92c) can be estimated as |^(p) — .^(p')| ~ where 
n ~ is the number of the levels available for the electrons. As a result, the 
largest term entering into the sums in Eq. (92c) is ~ r(*")T oc T^. It explains 
the power law dependence at low temperature. 

It is instructive to investigate the validity of the Wiedemann-Prantz law. Prom 
Eq. (99) we find for the Lorentz number 



L(r) ^ 3e^/^(r) 



/24 32 \ 5.%i 



— < 



192G^ I 



(100) 



The origin of the violation of the Wiedemann-Franz is the following. For the 
conductivity all the tunnelling pairs with the energies ~ T are roughly of 
the same importance. For the thermal conduction, however, the contributions 
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Fig. 10. (a) Lower and upper curves represent the conductivity a in the units of 
the Drude conductivity Goo for the Id and 2d cases, respectively [i.e., and 
exp((lnBo-))] as a function of ^j2dr{T) - temperature in the units of ^J6^T^l, as 
given by Eqs. (93), (94) in the fast energy relaxation appproximation (this curve 
does not depend on any parameters), (b) The relative Lorentz number L/Lq versus 
temperature in the units of ■sJZ^T^i for Id (upper curve) and 2d (lower curve) cases. 
Vertical dashed lines indicate the limit of the applicablity of the theory to the 
metallic region, see Eq. (80). 

of the low energy part is less important, and the best tunnelling pairs are 
different from those for the conductivity. It leads to the renormalization of the 
Lorentz number (100) by the factor of the order of unity. 



We conclude this section by noticing that the contribution of the rare pin- 
holes leads to the deviation from the "natural" assumption a oc C,ioc^^^'^\ 
cf. Refs. [25,28]. Wc think, however, that the relative contribution of such 
configurations at < T < T^i is a model dependent question, and do not 
pursue this line further. 



6 Stability of the insulating phase. 



Even though the condition (80) gives the sufficient condition for the metallic 
state to be stable, the consideration of the previous section does not prove 
the existence or the stability of the insulating phase, but rather gives the 
indications of the breakdown of the calculational scheme at T < Tj„ . To show 
the existence of the transition, we have to prove that perturbation theory from 
insulating side is also convergent under certain conditions. This analysis is the 
subject of the present section. 

The notion of the statisitical averaging is somewhat non-standard in the insu- 
lating phase and its general aspects are discussed in Sec. 6.1. To demonstrate 
the stability of the insulating phase we linearize the nonlinear equation (69a) 
in Sec. 6.2 and obtain an equation with the random coefficients somewhat 
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similar to that of Ref. [17]. Further consideration will be based on the sta- 
tistical analysis of this equation, using the technique described in Sec. 6.3. 
To begin with, we will investigate the simpler but yet instructive case of the 
zero-dimensional system 7 = to illustrate the general structure and make the 
connections with Ref. [16], see Sec. 6.4. To complete the study, we generahze 
the consideration to higher dimensions in Sec. 6.5 and calculate the transition 
temperature. 



6. 1 Averaging procedure 

Equations (69) form an infinite set of coupled non-linear equations whose coef- 
ficients are functions of the random energies and occupation numbers ni{e). 
Moreover, even though the equilibrium distribution n;(e) = tanh^^r nullifies 
the collision integral, the time of the relaxation of an arbitrary distribution 
ni{e) to the equilibrium one becomes infinite for the insulating phase. There- 
fore, ni{e) = tanh ^ has the meaning of the occupation number averaged over 
an infinitely long period of time (taken to infinity prior to setting the coupling 
to the bath equal to zero) and bears no information about the state of the 
system at a given instant of time. 

The basis in the space of many-body states of the system is formed by states 
represented as Slater determinants built on single-particle states (/,/?). Such 
Slater determinants correspond to the occupation numbers ni{e) = ±1. Thus, 
the decay rates F/ can be considered for each given set {rii}. Because the 
transition may occur only if the number of terms contributing to F/ is large, 
we can also perform the statistical average with respect to those occupation 
numbers. We assume them to be arbitrary with the only constraint being to 
fix the global energy of the system (microcanonical ensemble). Because the 
number of excitations involved into the formation of the decay rate is much 
smaller than the total number of the excitations in the system (which scales 
proportionally to the volume), averaging over the microcanonical ensemble 
with the energy E (counted from the ground state) is equivalent to the averging 
over the canonical ensemble with temperature T such that 

E - yi^VT^, (101) 

where u is the averaged density of states per unit volume and V is the volume 
of the system. We assume T » MS^ and verify later that the transition occurs 
only at such temperature. 

To find the distribution functions of the local and thus strongly fiuctuating 
quantities, we perform the ensemble averaging on the final stage of the calcu- 
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lation. Thus, the averaging procedure in this section is defined as 



NS(/2 



hm 



/ nn^ E 



2 cosh 



2T 



(102) 



where the number of levels N and the domain of integration is in accordance 
with Eq. (48a). Formula (102) docs not take into account the level repulsion, 
see Eq. (48b). This repulsion will be included later when it is necessary. 

Let us notice that the same averaging formula would arise if one considered 
the probability of the occupation (and not the average occupation) of the 
levels formed as a result of the arbitrarily weak interaction with equilibrium 
phonons kept at temeperature T. In this respect temperature T in Eq. (102) 
has a meaning of the experimentally measurable quantity. 

Having discussed the issue of how to average, we are prepared to calculate the 
characteristic function (66) to establish the conditions of the stability of the 
insulating phase. 



6.2 Linearized ImSCBA equations 



To begin the actual analysis of the insulating phase, we notice that in the 
absence of the external bath vf"'^^\e)^ substitution r/(e) = solves the self- 
consistency equation (69a). However, in accordance with the discussion of 
Sec. 4.2 the order of limits is important, so that we have to investigate the 
stability of the solution r;(e) = with respect to small but finite coupling 
to the bath. The smallness of r;(e) oc h enables us to linearize the spectral 
density in Eq. (69a) as 

Me)^5{e-ii) + l-^^^^, (103) 
^ (e - 6) 

where the second term is understood as the principal value 

We substitute Eq. (103) into the equation for r;(e), see Eq. (69a), notice that 
the probablity of matching the levels into exact resonance equals to zero, and 
keep the terms only linear in r;(e). Taking into account that, according to 
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Sec. 6.1, ni{e) — ni — ±1, we obtain 



a)] 



(i,o 



+ E 



:i05a) 



X 



2r^,(e-e/,+e/3) + r/3(e/i+e/2-e) 



where the denominators are defined in the sense of Eq. (104), and 



li,l2;h 



1; if ni^ = = -71^3 



±1; 



0; 



otherwise. 



The decay due to the connection with the bath 



^(bath) 



coth 



2T 



+ ni 



(105b) 



(105c) 



is a smooth function of e. As before, the coordinate p is assumed to be the 
same in all terms in the equations unless it is specified explicitly otherwise. 

To deal with certain superfiuous logarithmic divergences in the future analysis 
of Eq. (105a), we introduced an ultraviolet cut-off with the help of the function 



0a{x) 



1; \x\ < A 
0; \x\ > A. 



(106) 



This ultraviolet cut-off is introduced for technical convenience only and 5^ -C 
A < MSi^, compare with Eq. (86). The contribution from the part excluded by 
cut-off does not contain small denominators. This part is self-averaging, pro- 
portional to b{uj), and, thus, it is not relevant for the question of the stability 
of the insulator, see Eq. (65). 

The second term in Eq. (105a) is the effect of the one particle tunneling into 
the neighboring localization cells. Due to the condition (48f) this term alone 
does not lead to any significant effects. As this term corresponds to the solution 
of the one-particle Schrodinger equation, it does not depend on the occupation 
numbers rii. On the other hand, the last term in Eq. (105a) describes the decay 
due to the excitation of electron-hole pairs. The availabihty of orbitals for such 
a process is controlled by the set of ri/ = ±1. 

Three following subsections are devoted to the statistical analysis of Eq. (105a) 
with respect to the distribution (102). Namely, we formally solve Eq. (105a) 



52 



as an expansion in small parameters A^, P: 

and calculate the characteristic functions 

lyM(s) = ^cxp (-srf ''"^(e))) . (108) 

for each term. The insulator is stable, see condition (66), if the typical value 
of each term decreases, i.e 

lim lim lim W^''^"'\s) = 1, (109) 

6— >0 V— >oo fc,m— »oo 

for any fixed s > 0. 

6.3 Mayer-Mayer cluster expansion for the characteristic functions 



In the actual calculation of the characteristic function (108) we will use the 
analogy of this function and averaging procedure (102) with the partition 
function of the system of classical particles interacting through a certain many- 
particle potential. The energy plays the role of the coordinate, occupation 
number = ±1 is equivalent to the particle spin, and parameter s is analogous 
to the "inverse temperature" . In order to find these interaction potentials we 
introduce the diagrammatic representation of Eq. (105a). This representation 
and all of the notation are shown on Fig. 11. 

It is obvious that k + m iterations of the diagrammatic equation shown on 
Fig. 11 produce diagrams for each term r*^'^'™^ of Eq. (107). Every term depends 
on all the orbital energies and all the occupation numbers. At the same time, 
it can be represented as a sum of terms depending only on the set {^i, ni}'l^'l^"^, 
such that 

r(^,-) ^ Y: U(^i„ni,; . . . ; 63.+., m^.^J, (110) 

where we omitted the spatial coordinate p. Thus, averaging (102) in the ex- 
pression (108) is indeed equivalent to the calculation of the partition function 
of iV — > cxD particles interacting via (3A; -|- m)-particle potential U. 

To perform this calculation we employ the procedure known as Mayer-Mayer 
cluster expansion [29]: 

00 

lnp^(fc,m)^ ^ lnVFp(^'™), (111) 

p=3k+m 
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e + Ci3 = ei+ e + ea = 61 + 62 

Fig. 11. Diagrammatic representation of Eq. (105a). Solid lines bear an orbital 
energy ^; and an occupation number ni = ±1. Double lines bear an orbital energy 
^i, an occupation number n/ = ±1, and an energy e. The energy carried by the 
double line is conserved in each hopping vertex (bold rectangle) and changes in the 
interaction vertex (white square). The algebraic sum of the energies entering the 
interaction vertex through the double lines and the orbital energies entering through 
the solid lines is conserved (signs arc determined by arrow directions). All the other 
rules are defined on the figure, and 6^ is defined in Eq. (106). 

where the summation is performed over linked clusters of p particles. Analo- 
gously to the calculation of the partition function of a classical gas interacting 
via a two-particle potential U{ri — r2), whose leading term of the cluster 
expansion is given by 

e T —1 ^ 

for each diagram (110) we introduce the cluster function / = e~*^— 1, as 
shown in Fig. 12. On the same figure we show the diagrammatic representation 
of the averaging procedure and linking the clusters. 

This procedure is especially suitable for the present problem as it systemati- 
cally takes into account contributions from regions of the phase space where 
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lnW2{T) 



d r 



(a) 



_ ^-8NuU{ii^,ni^;...;ii^,nij,) _ 



(b) 



(c) 



x/s(---;6,^z;---) 



\ 



d6 e 2T 



/ ^ n,=±l ^ ^"^^^ 2T 

x/(...;6,n,...) 

x/^(...;e«,ni;...) 

Fig. 12. (a) Definition of the cluster function /; Njj is the degeneracy of the diagram 
- the number of fine permutations keeping the diagram intact, (b) Independent 
averaging over the position and the occupation ni of a level I. (c) Connected 
average of two different cluster functions / and /' (an arbitrary number of cluster 
functions can be averaged analogously). The contribution of each linked cluster to 
InWj^*^'™^ (p being the number of levels averaged over) should be divided by the 
number of symmetries of the averaged diagram. This coefficient will be written 
explicitly in each specific case. 

the interaction potentials U assume large values, «.e., of small resonant denom- 
inators responsible for delocalization. The next two subsections are dedicated 
to the application of this procedure. 



6.4 Statistical analysis for zero-dimensional system. 



In this section we consider only transitions inside one localization cell, m = 0, 
or, in other words, 7 — > 0. It suffices for our purposes to limit ourselves with 

To understand the most crucial features of the cluster expansion, we will 
consider a few lowest order terms explicitly, and then analyze an arbitrary 
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order term. It will be useful for us to introduce the notation 

C::C'' = E4-Ee.. (113) 

i=i j=i 

for the energy of an excitation consisting of Ne electrons and Nh holes. 
The lowest order terms are 

r(O'O) = T6(0); (114a) 



(114b) 



where 6'a is given by Eq. (106). In Eq. (114b), we chose not to join the electron 
decay and the hole decay even though they are equal to each other. This is done 
to keep the structure of the subsequent terms more transparent. Moreover, 
this form will be useful for the discussion of the modification of ImSCBA in 
Sec. 7.3. 

Let us use the machinery of Sec. 6.3 to find the characteristic function of the 
rate (114b). The corresponding three-particle potentials are given by 

C/l%3 = '2U^2;3 = '^3,2 . (US) 

switched on for certain combinations {t^i, n2, ^3}, see Eq. (105b). The cutoff 
9a{€ — 5^2) will be incorporated in the definition of cluster function / below. 

The main conribution comes from clusters with p — 3 particles, see Fig. 13a: 



(1,0) 




d6 



k n^il2cOsh|. 



E 



/l2;3; 



/l2;3 — /l2;3 + 2'^12;3 + 2'^12;3 {/21;3 + /l2;3 ^ + /21;3 } 

/r2;3 = (e-^^-^^ -1) OA{e - ^l,) , 
/i'^2;3 - (e-""-- -1) e^{e - SL) , 



(116a) 

(116b) 
(116c) 
(116d) 



We assume that b{uj = 0) > 0. This assumption does not correspond to any 
physical phonons. Actual frequency dependence b{uj) oc lu^ is important if one tries 
to calculate the physical conductivity in the insulating regime. It is however not 
important for the determination of the stability of the inslulating phase. 
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Fig. 13. The 3-, 4-, and 5-particle contributions to Inl^^-'^''^). For three particles 
(a), the cross-correlations between the electron and the hole contributions, the last 
term in Eq. (116b), are shown in the second and third lines. For 4 and 5 particles 
(b,c) only the hole-hole terms are shown. When not displayed, the degeneracy for 
all diagrams Nu are the same as on panel (a). 
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where the coefficient 1/2 and extra factor of Njj = 2 in the exponential in 
fi2:i take care of the symmetry /{2.3 = /2i-3- second term on Fig. 13a. The 
last termm in Eq. (116b) describes the cross-correlations between the electron 
and the hole contributions as the potentials [^"^(Ci, C2, Cs); ^^(G, {1, Cs); and 
U^iii, 6, 6) diverge for the same sets of 6, 6, 6- 

Performing the integration with the help of Eq. (105b) and definition of Y , 
see Eq. (69a), we obtain using the condition (112): 



InW 



(1,0) 



3 



where the numerical coefficient l/\/3 appears due to the cross- correlation 
terms in Eq. (116b). 

Parameter TM.^ 2!^^ roughly corresponds to the total energy of all electron- 
hole pairs participating in the process of the decay of a quasiparticle with the 
energy e. The explicit expression for A^d is 

M{y,z) = y[ ^"^^"^ - T^(yx.)| cosh. 

2 J cosh(xi +X2 + z) ni=i,2 cosh(Xi + z) 

Performing integration in Eq. (118) with the help of Eq. (48h), we obtain for 
T > M6^ 

M = 2M, (119) 

independently of T and e. This occurs because of the additional restriction on 
the phase volume for the decay of one-particle excitation into the three-particle 
excitations by the energy dependence of the matrix elements. 

Although the metal-insulator transition is a property of large orders of pertur- 
bation theory, see Eq. (109), its precursor can be seen already in the character- 
istic function of the lowest order term (117). Indeed, condition InW'^^'^^s*) ~ 
1 gives the most probable value P^^^^ ~ 1/s* of the distribution. Requiring 
rllf < r(°'°), we obtain the condition T < T*, where 

fXMT, 



~ 1, (120) 



In fact, using Eq. (115), one can rewrite 

/12;3 = ^(1 + /f2;3)(l + /ll;3)(l + /l'2;3) = ^ [exp(-6C/'^) - 1 

i.e., join all the term in one cluster function. This joining, however, would obscure 
the structure of the higher order terms and we choose not to perform this operation. 
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i.e., T* is roughly of the order of the temperature Tj„ hmiting the stabihty 
of the metalhc phase, see Eq. (80). Consideration of higher orders of the per- 
turbation theory to be performed shortly will only slightly refine expression 
(120), see Eq. (131). 



To complete the analysis of the lowest order term (114), we verify the va- 
lidity of the cluster expansion, which amounts to finding the condition for 
which \ lnW^'^'^\s)\ > \ lnWiY\s)\. Diagrams for lnWiY\s) are shown in 



Fig. 13b, c. They give 



X 

4 



7 (/l2;3/l2;4 + fl3;2 fu;2 + 2/l3;2/24;l + 4/i2;3/l4;2) , 
^ ^ ^ (121) 



2T 



,.1 V .,tli2coshi 



X ^ (/l2;5/34;5 + 4/5l;2/53;4 + 4/i2;5/53;4) , 



Using /i2,3 from Eq. (116b), we can estimate these expressions as 



where M is defined in Eq. (48h). Comparing Eq. (122) with Eq. (117), we find 



lnPFi^'°) 



> 1. The 



that for T ^ the cluster expansion is justified even for 
origin of this suppression of the larger clusters is the same as the one controlling 
the virial expansion for the classical gases - clusters involving more particles 
impose more restrictions on the phase volume. 



Our next task is to consider an arbitrary order term r'^'^'^-'. To explore how the 
lowest-order result Eq. (117) is modified, we consider F^^'^) first. Performing 
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Fig. 14. The leading contributions to the cluster expansion 

of lnW^(2,o)_ The cross- 
correlations between the electron and the hole contributions, second and third lines 
in Eq. (125b), are not shown. 



one more iteration in Eq. (105a) with the help of Eq. (114), we find 



p(2,0) _ p(0,0) ■y4r4 V V 



each term is shown on Fig. 14. Here prime in the second sum excludes the terms 
with /4, = li, I2 and /g = ^3, see Fig. 7b. 2. (Let us note that the placement of 
crosses on Fig. 14 automatically takes this exclusion into account.) Similarly 
to the expression for F'^^'^'' in Eq. (114), we chose not to join some similar 
terms together. 
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The corresponding 6-particle potentials are given by 



eh 



hh 



(e - El,f (6 - En,f 



(124) 



The leading term of the cluster expansion shown on Fig. 14, is given by [cf. 
Eq. (15)] 



/ A2;3\ _ fee A2;3\ , ^ reh /12;3\ , ^ r/ie /12;3\ , ^ /-/i/i /12;3\ 

I 1 fee /12;3\ / /ee /12;"^ 
"•"2'' V45;6j \/ V54; 

1 



3\ I reh /l2;3\ 

1 f/ie / 12; 3\ f |-/ie / 12; 3\ , rhh ( 12; 
4/ 1,45; 6;\/ 1,54; 6;+/ 

V ,,„„/l2:3\ ^ 
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^a(6 - 5 


L) - 


/•/le /12;3\ 
/ 1,54; 6 ; 




"^-2.t/'^=(i2;3^ 


-1 


^a(6 - 5 


L) - 


fhh /12;3\ 
/ 1,54; 




_4 rrhh/'l2;3> 

e V54;6y 


-1 


^A(e - E 





2 + r(g^)]} 



(125a) 



(125b) 



(125c) 



Similarly to Eq. (116b), the second and the third lines in Eq. (125b) express 
the cross-correlation between different terms in the lowest order perturbation 
theory. These correlations do not proliferate into higher order terms; in par- 
ticular, U^^ from Eq. (124) is not invariant under permutation 1 <-> 2. 

Substituting Eq. (124) into Eq. (125) and using the approximate expression 
for the integral 



2-„{y, z) = [ dxi... dxn 
Jo 

lnz"l/-V2 + C'(i) 



exp 1^ 

m-i 



y 



x\...xl^ 



- 1 



(126) 



(n-l)\ 



-; y < 



2n 
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we find 



lnW^(2,o)^__L(^5r(0'0)) 



1/2 



'6XMT' 


I 

X In 


[ k \ 





(127) 



where A is the scale introduced in Eq. (106), and the coefficient l/-\/3 comes 
from the cross-correlations in Eq. (125b). 

The analysis of the corrections to Eq. (127) is performed similarly to the 
estimate of Eq. (121), see Appendix C. It shows that their main effect is the 
more accurate determination of the cut-off of the logarithm such that 



So 



(128) 



so that Eq. (127) is valid in a leading logarithmic approximation. This value 
of the cut-off is in agreement with Ref . [20] . 

All the higher order terms are considered analogously, and we find 

A" 



/7rsr(0'0) f6\MT\ 



3 V y (^-1)! 



In 



n-1 



(129) 



To investigate metal-insulator transition we have to consider the behavior at 
n — > oo. We obtain from Eq. (129) with the logarthimic accuracy: 



1 -1 



CJ 
(130a) 



Equations (130) constitute the central result of this section as they describe 
the complete statistics of each term in the ImSCBA series. If many particles are 
in the excited state (T is finite), formula (130a) predicts the phase transition. 
Indeed, the stabihty criterion (109) is violated at T > T*, where 

^-^ln- = l, (131) 

where we used the replacement (128) in the argument of the logarithm. Equa- 
tion (131) refines Eq. (120). 

This would mean that insulator becomes unstable at T > T*. This conclusion, 
however, is an artefact of the averaging procedure (102) which neglected the 
one-particle level repulsion, see Eq. (48b). To take this level repulsion into 
account we use the following qualitative consideration. Inspection of each term 



62 



in the pcrturbaiton theory, see, e.g., Eq. (15) or (125), shows that W^^'^^ is 
contributed by 2n integrations over — T < < T. It means that is 
essentially determined by the probability to find 2n levels within the interval 
with the width T. According to Eq. (48b), this probabihty is suppressed by 
the level repulsion. As a result, Eq. (130a) is modified as 



l^W^n»m ^ _ (,r(o>o))-/^ I In AV exp 




, (132) 



where function V is defined in Eq. (48b) , and c is the numerical factor of the 
order of unity. Because \imx_^^x~^V{x) — oo, we find 

lim lnW^("'°) = for any T. 



Equation (132) describes the impossibility of a genuine phase tranistion in a 
finite system where the unphysical packing of infinite number of levels in a 
finite enery strip is forbidden. 

Even though for a zero-dimensional system temperature T^, does not have a 
physical meaning of the transition temperature, the finite coupling between 
the locahzation cells leads to the phase transition with Tc very close to T*, as 
we will show in the next subsection. 



6.5 Statistical analysis and metal-insulator transition in finite- dimensional 
systems. 



The goal of this subsection is to generalize Eqs. (130) by including hopping 
into the neighbouring localization cells. As a result, the effect of the finite size 
(132) will be overcome and the genuine metal-insulator tranistion will occur. 
Similarly to the previous subsection, we will consider explicitly a few lower 
orders of the perturbation theory, and then obtain the result for an arbitrary 
order. 

The contribution which involves hopping only (its first two terms F^'''^^ and 
r'^*''^^ are shown in Fig. 15) corresponds to the purely single-particle problem. 
Since the exact single-particle eigenstates are assumed to be localized, see 
condition (48f), this contribution does not lead to delocalization. Thus, we 
start with the term r^^'^\ which involves tunneling of the particle and creation 
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of an electron-hole pair. By iterating Eq. (105a) twice, we find (see Fig. 15): 



h,...,l4;a 



2 / n 

e - E? 



X 



+ 



^'3 

hit 



^ -hit 



+ 



^hL^h,h;h^^{^ ~ -hh) ^^[-hh - ^ 



(133) 



where we used the short-hand notation (113), and = ii{p), Ci+ = 6(P + 
Prom Fig. 15 we find (omitting overall numerical coefficient) 



-) 



(134) 



Comparing Eq. (134) with Eq. (127), we see that by replacing the electron-hole 
pair creation to a single-particle hopping, we always introduce the additional 
smallness. Thus, we can anticipate that the number of hoppings must be as 
small as possible to maximize the overall W^"^''"'^ for fixed m + n. The only 
reason to include hopping at all is to overcome the finite number of levels 
restriction which suppressed the transition in 0-dimensional case of previous 
subsection. 

The result for the arbitrary order perturbation theory may be immediately 
obtained by examining the differences between Eq. (134) with Eq. (127), and 
analogy with Eq. (129): 



0) 



1/2 (qxmtY 



X 



C(n, m) 
(n -|- m — 1)! 



In^ 



n+m—l 



(sr(O'O)) 



-1/2 



(135) 
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InW, 



(0,2) 



(1,1) 




Fig. 15. The leading contributions to lnW^^^^\ InVF^^-^), and lnl^(^'^). The 
cross-correlations between the electron and hole contributions are not shown. 

The most important difference with Eq. (129) is the presence of the additional 
factor C{n,m) which has the meaning of the number of ways to order the n 
interactions and m hoppings with respect to each other. The expression for 
these coefficients is 



m CO 



C(n,m) = Jl Sn,j:N,cxp 

i=l Ui Ni=0 i 



T 



V 



, (136) 



where pg is an arbitrary site on the lattice, a are the vectors connecting each 
site to its nearest neighbours. The meaning of the coefficients, C{n,m), is the 
total number of random lines consisting of m segments on the d-dimensional 
cubic lattice. Integer Ni is the number of electron hole pairs emitted by the 
electron between the hopping event i and the hopping event {i + 1). We in- 
cluded the effect of the level repulsion within each localization cell in the same 
spirit as in Eq. (132). The argument of V counts the total number of electron- 
hole pairs within localization site p and takes into account the fact that the 
random path may traverse one site more than once. 

To investigate the stability of the insulating state, see Eq. (109), we have to 
study the limit of m = m*(n), n ^ oo, where m*{n) determines the direction 
to maximize InVF^"'"*^ for fixed n. 



As before, we will restrict ourselves by the condition 

7 ~ A < AM. 



(137) 
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Simple counting of combinatorial factors in Eq. (136) shows that under con- 
dition (137) maximum is achieved for the maximal possible number of the 
electron-hole pairs per a localization cell, iVj ~ T/5i^. It translates into esti- 
mate 

m*(n) ~ (138) 

Substituting Eq. (138) into Eq. (136), we find C[n,m*(n)] < (2d)'"*(") and 
obtain from Eq. (135): 

where Ci{d) is the number of the order of unity which we were not able to 
calculate, and all the other entries are the same as in Eq. (130a). 

Applying stability criterion (109), we arrive to the main conclusion of this sec- 
tion - insulating state is stable only for T < Tc, where the critical temperature 
is given by 

As / is not exponentially small, / ^ e~^^^^ , Eq. (140) may be expanded as 

n^T,[l + 0{XM\lnI\)], (141) 
where is given by Eq. (131). 

At temperatures larger than Tc metallic phase is formed. Together with the 
material of Sec. 5 proving the stability of the metallic phase at T > Tj„ ^ Tc, 
see Eq. (80), this completes the proof of the existence of the metal insulator 
transition. 

Kinetics of the system near the transition itself is a complicated problem which 
we hope to address in a separate publication. However, some conclusions can 
be drawn already from Eqs. (135), (136), and the estimate (138). Indeed, 
we concluded that the best paths are those that maximize the number of 
electron-hole pairs in a given localization cell before hopping to a neighbouring 
one. It means that the self-intersections of the random path in the Eq. (136) 
are forbidden and the spatial part of the problem becomes equivalent to the 
statistics of the self-avoiding random walk. 

The latter observation enables us to conjecture the critical behavior of the 
spatial localization length ^ at T — > — 0. The latter length is defined from 
Eq. (136), as 

i typ 
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where the configurations giving the largest statistical weight in Eq. (136) are 
meant by typical. The correlation length in the Fock space (typical distance 
between resonances) diverges in the vicinity of the transition as [22] 



^typ — 



Tc-T 



(142) 



hence, the number of segments m ~ ^^/jTc — r| ^. We thus conclude 



where is the correlation length index for the self-avoiding random walk in 
d dimensions, = 1, 1^2 = 3/4, 1/3 = 0.59 . . . , u^^^ = 1/2, see Ref. [30]. 



7 Validity of ImSCBA scheme 

This section is devoted to the analysis of the contributions of the processes 
not taken into account in our ImSCBA calculational scheme, see Figs. 4,5, 
and Eq. (68). The latter approximation corresponds to neglecting the level 
shifts due to the electron-electron interaction. The former will be shown to 
correspond to the renormalization of the electron-electron interaction by inter- 
mediate virtual processes, and to certain interference effects. We will consider 
the effects of those contributions separately in the following subsections. 

We note here, however, that each remaining diagram, denoted by U, is a 
random quantity. Therefore, we will have to analyse the statistical distribution 
of the remaining diagrams and check that their distribution functions have the 
scale parametrically smaller than the scale of r;(e) calculated in the previous 
section. The signs of the majority of the diagrams are random, so it is more 
convenient to use the characteristic function, see, e.g., Eq. (108), in a Fourier 
transform form 



where the averaging procedure is defined in Eq. (73). To calculate function 
(144), all the machinery of Sec. 6.3 is applicable after the replacement s — > —iq, 
and U on the Fig. 12 should be understood as an analytic expression for the 
real or imaginary part of the corresponding diagram. 

We will present in detail the analysis for the insulating phase in Sees. 7.1 - 
7.3; the corresponding consideration for the developed metallic phase is simple 
and is summarized in Sec. 7.4. 




(143) 



(g) = (exp (iqU)) 



(144) 
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7. 1 Effect of the interaction renormalization. 



To understand the origin of the interaction renormalization, let us consider 
the contributions c5,c6 from Fig. 4, see also Fig. 16. These contributions are 
the only third-order terms that may lead to the finite decay rate - all others 
are either insignificant corrections to the Hartree-Fock potential (c2,c3), or 
(cl,c4) the first order Hartree-Fock shift in the second order diagram (bl,b3). 
Direct calculation of the diagram (c5,c6) yields 



Im Ef (^) = ^ E / ■ ■ ■ d65 A, (eO . . . A, (e,) 



ll,...,l5 



^ih ^hk ^kh^^ (e - ei - 62 + €3) i'i,,i^.i^[ei, €2, €3) 1^- 

C2 ~ £3 — €4 -|- €5 

+ (2 ^4; 3 ^5) 



j.A(Fig4c6)(^^ ^ 1 ^ j d65 ^,(60 . . . A,{e,) 



h,...,h 

X 



+ (1^4; 2 ^5) 



£1 + £2 ~ Q ~ ^5 

(145) 



where the notation is the same in in Eqs. (69), and P denotes the principal 
value. It is important to notice that using ImSCBA Green functions in this 
expression rather than the bare Green functions is not an overstepping of the 
accuracy of the calculation. 

Expression (145) illustrates the well-known principle of constructing higher 
order contributions from the lower ones. Namely, to obtain the imaginary part 
of any contribution to the self-energy one can cut the diagram in all possible 
ways [cuts are shown by vertical dotted lines on Fig. 16]. The cross-section 
produces a 5-function for the energies of the particles crossing the cut. The 
two parts of the diagrams on each side of the cut correspond to transition 
amplitudes. Being real, they can be incorporated into the redefinition of the 
constants of the initial Hamiltonian (for the clean Fermi liquid it was first 
realized by Eliashberg [31]). 



For example, expression (145) can be obtained from Fj" (e) of Eq. (69a) by 
the replacement of the bare interaction matrix element Vi^i^ with the poten- 
tial renormalized by excitation of virtual particle- hole pairs, SVeh-, or particle- 
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particle pairs, 5Vee, see Fig. 16b: 



[^Ke]!i = ^E / dea dee^,(.5)^.(ea)^it^a-P :f f ^ "'^^'f ; 

4 J ei + 62 - £6 - ea 

where we are using the short-hand notation li = {li,ei) in the first fine of the 
expression. 

For the insulating regime we use the leading term in Eq. (103), Ai{e) ^ 
S{e — ^i), and calculate the characteristic functions W^^'^'^q) and W^^^'^{q) 
as shown on Fig. 16c. Assuming |e2 — Cs], |e2 + ^il ^ 5^, T > M6(^, we find 



-2n\q\X^T; |e2 - £3! < MS^; 

(147a) 

-; \e2-e^\>M6^- 



n\q\X'^M5^_ 



In Vl^'^^-(g) = i. (147b) 



Equations (147) describe random quantities with the characteristic scale of 
the distribution I^V^It^^p given by the coefficient multiplying This width 
should be compared with the bare value of the interaction constant A5^. For 
T <T^, [see Eq. (131)], we find 



\SV\,,, ^ 1 / T 



XS M \T* 



(148) 



i.e., the typical value of the correction to the interaction constant is paramet- 
rically smaller than the bare value. On the other hand, the distribution of the 
SV has the same algebraic decay at large values as the distributions of the 
ImSCBA quantities F, see, e.g., Eq. (129). It means that substitution of the 
renormalized constant SV in, say, second order perturbation theory formula 
(114) will produce the distribution function of the form as the fourth order 
perturbation theory result (125), but with the coefficient smaller at least by 
the factor of M. The same is true for any order, and therefore, the interaction 
renormalization (147) can produce only perturbative in 1/M corrections to 
the value of the transition temperature T^.. 

It is not difficult to see that certain cross-sections of higher-order diagrams 
may be ascribed to the higher-order corrections to the interaction vertex, see 
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b) 





In W^^-> 



-iq 




-n 



Fig. 16. Reduction of the third order diagrams (a) to the renormahzation of the inter- 
action constant in the SCBA scheme (b). Calculation of the characteristic function 
(c) using notation of Fig. 12. 

Fig. 17. However, not all of the cross-sections can be taken into account in a 
such a fashion, see, e.g., cross-section (c2) of Fig. 5, or Fig. 18a. Remaining 
terms describe the effects of the particle permutations in the final state which 
will be discussed in the following subsection, see Fig. 19, and generation of the 
interaction vertices involving a larger number of the particles, see Fig. 18c. 



Statistical analysis of the higher-order corrections to the vertices is performed 
in the same fashion and produce distributions similar to that of Eqs. (147) with 
the smallness (148) in higher and higher powers. We thus conclude that the 
vertex renormahzation does not lead to any dramatic effects; it only produces 
perturbative corrections to the transition temperature Tc calculated within 
ImSCBA scheme. 



7.2 Effect of the particle permutations in the final state. 



Cross-sections of the fourth-order diagrams, not included in Fig. 17, are shown 
in Fig. 19 (5.1a-5.2d) together with the two ImSCBA cross-sections (5.1 and 
5.2). One can notice that all of these cross-sections contain the sum over the 
same final 5-particle states (the same orbital indices of the particles and holes 
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2 ^1- 



+ 0{X^ 



+ 

y ♦ 







Fig. 17. Reduction of the certain cross-section of fourth order diagrams, Fig. 5, to 
the renormahzation of the interaction constant in the SCBA scheme. 

crossing the cut). Moreover, the analytic expressions corresponding to the two 
parts of a cut ImSCBA diagram for a given final state (transition matrix el- 
ements) are equal. For non-ImSCBA diagrams the two resulting expressions 
arc different, but they always can be reduced to the transition matrix ele- 
ments for the diagrams 5.1 and 5.2 by a permutation of the orbital indices 
of the particles in the final state. This observation enables us to identify the 
non-ImSCBA cross-sections 5.1a-5.2d as interference terms in the transition 
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(d) 





□ 








□ 



-iq 



Fig. 18. a) Reduction of a crosssection of the 6-th order diagrams to the three particle 
interaction, b) SCBA diagram generating the same state in the crosssection. c-d) 
Comparing of the distribution functions of the non-equivalent blocks in the SCBA 
diagram (d) and in the 3-particle interaction diagram (c). 



probability, while ImSCBA corresponds to the replacement of the square of 
the sum in Fig. 19 by the sum of squares of individual terms. 

This consideration for the fourth-order self-energy illustrates the general rule 
of construction of all lowest-order diagrams containing a cross-section of elec- 
trons and Nfi holes. Let us assume for a moment that electrons and holes are 
distinguishable particles and introduce the quantum-mechanical amplitude 

Ai{{li}'iti', {{'m-i}i^=i) of the transition to the given final state for a given path i 
(represented by a half-diagram with a definite assignment of orbital indices to 
the free electron and hole lines). The decay rate can then be represented as 



i=l 



ll 



1 



1=1 



i Pl,Pm 



(149) 



The factorial prefactor takes into account the fact that different permutations 
of the indices {li},{mi} correspond to the same final state, the sum over 
the permutations represents the usual antisymmetrization of the quantum- 
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Fig. 19. Cross-sections of the fourth-order diagrams for self-energy, not included in 
Fig. 17. The orbital indices of the particles crossing the cut {Ii,l2,l3,mi,m2) are 
summed over. For fixed values of these indices each of the two parts of the diagram 
with the free ends removed corresponds to the partial transition matrix element. 
Summation of the latter over all permutations of the orbital indices of identical 
particles in the final state gives the total transition matrix element for a given final 
state. The permutations Pi and Pm arc independent; they act on electron {li,l2,h) 
and hole (mi, 777,2) indices, respectively. The decay rate is then given by the square of 
the total matrix element, summed over all final states. Since different permutations 
of electron and hole indices give the same final state, the factor l/(3!2!) in front of 
the sum is necessary. Diagrams (2) and (4) of Fig. 7c are also generated for / I' 
or m" ^ m'. 
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{l2,l';mi) 




{Ii,l2,h]rni,m2) 



{li,l2\m') 



Fig. 20. Different paths for tlie transition from a one-particle state I to a five-particle 
state (li, 12, h', mi, 1712), corresponding to different intermediate three-particle states. 
The state {l2,h,h','fni,'m2) is identical to {Ii,l2,h','mi,m2)- 

mechanical amplitude . 

Similarly to the consideration of the quantum interference effects for a single 
particle in a disorderd potential, the resulting double sum in Eq. (149) can be 
separated into diagonal and off-diagonal part 



The first term in Eq. (150) is nothing but the ImSCBA series corresponding 
to the summation of the probablities of the paths. Random-sign last term is 
the quantum interference contribution to the particle lifetime. 

Even though the number of the interference terms is much larger than the 
number of diagonal terms one can still argue that they do not affect the value 
of the transition temperature. 

Indeed, if the number of the relevant terms were large, one would be able 
to apply the central limit theorem for the quantum-mechanical amplitudes 
rather than for the probabilities. It would result in a distribution function of 
the same scale as the one obtained by the diagonal approximation but of a 
different shape {e.g., for the large number of statiscally independent ampli- 
tudes Porter-Thomas distribution would replace the Gaussian one). However, 
the characteristic functions obtained in Sec. 6.4 arc non-Gaussian (as they 
are not analytic at s — > 0) and have long algebraic tails at large values of F. 
This indicates that the result is contributed by the largest term in the sum. 
Therefore, the interference contribution affects the distribution function in the 
range of the most probable values but not the tail of the distribution. As the 



It is easy to see that the diagrams (2) and (4) of Fig. 7c are also included in 
Eq. (149) 



n oc ^ \Ae{{li}, {m,})f + (off-diag.) 



(150) 
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transition temperature is controlled by the tail, the interference term is not 
important for the position of this temperature. 

To quantify the qualitative consideration above we evaluated the characteristic 
function for the quantity on the right-hand side of Eq. (149), calculated in the 
nth order of the perturbation theory (n ^ 1), 

ln/exp|-s5^2 J2 7K{^i}> Wi}) I ) = -v^ (1 - V^amt + ■ ■ ■ ) ; 



{li},{mi} 

2 



(151) 



Derivation of Eq. (151) is relegated to Appendix D, and we neglected In A 
factors in the expression for aint- 

Leading at s — > term is none but the the diagonal contribution calculated in 
SCBA approximation. The subleading term is the interference term, and the 
parameter a^nt has a meaning of the relative contribution of the off-diagonal 
terms, affecting the distribution at small values of 7. Values of aint contain the 
n! factor, so it apparently becomes large no matter how small the interaction 
constant A is. However, as we discussed in the derivation of Eq. (132), the 
order of the perturbation theory involving only one localization cell is limited 
from above, n < n* ~ T/6(^. As the result for 7 = 0, n = n* we estimate 



XT] ^"-^ 



K J 



M-^"' < 1, (152) 



which means that the interference processes do not affect the most important 
part of the distribution functions even in the largest possible order of the 
perturbation theory without hopping between localization cells. 



One could think that the inclusion of the hopping would allow the growth of 
the aint beyond the estimate (152) at n > n*. However, it is not the case. 
Inclusion of tunneling into neighboring localization cells will lead to further 
suppression of the interference effects. The presence of tunneling vertices on 
the diagram makes it impossible to interchange particles residing in different 
localization cells, compare Figs. 19 and 21. It is possible to check that keeping 
the combinatorial n! in the perturbation theory involving n particles in the 
final states spread over m <^n localization cells would require the tunneling of 

Moreover, experience gained in the study of the critical behavior of the Anderson 
transition on the Cayley tree [22] suggests that the transition itself is associated 
with the reconstruction of the tail of the distribution. Therefore, we do not see any 
reason to believe the the interference terms can affect the critical behavior in our 
problem either. 
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Fig. 21. Suppression of the exchange processes shown on Fig. 19 due to the inclusion 
of the tunnehng vertices. Diagram a) has only one exchange counterpart b) whereas 
the diagram c) does not allow permutations at all [cf. diagram (5.2) which produces 
three exchange counterparts]. 

all n particles by the distance of the order of m. Each tunneling event brings 
additional smallness /, and thus we estimate 

where we used m ~ m* from Eq. (138). Thus, no accumulation of the factorial 
terms is possible, even in the vicinity of the metal-insulator transition, and 
the SCBA calculational scheme is valid not only on qualitative but also on a 
quantitative level. 



7.3 Effect of the single-particle spectrum renormalization. 

We now turn to the study of the effects on the interaction and tunneling which 
may be viewed as a change in the properties of single-particle excitations. This 
change includes the level shifts S^i and the possibihty of mixing with other 
orbitals (which is present even without interaction due to the hopping /). 
The level shifts are encoded in the statistics of the real parts of the diagonal 
part of the self-energy ReSf which were neglected in ImSCBA approximation, 
whereas the bubbles non-diagonal with respect to the initial and final states, 
see, e.g., diagram (c3) of Fig. 7, are responsible for the mixing. 



nAexp 



n\ In I\6(; 
T 



(153) 
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The lowest order diagram is the Hartree-Fock potential of Fig. 4a 

^wf^-^EK-sg^e^al^S, (154) 

where we subtracted the value of the potential at T = 0. The latter term is 
assumed to be incorporated into the one-particle Hamiltonian. The character- 
istic function (144) is found using the averaging procedure (73) and Eq. (48g) 
to be Gaussian 

Iniy^"^ = -e\\T\^2 ~ -qHl (^A^ {L^ . (155) 

This means that the variation of the Hartree-Fock potential with the distri- 
bution function is much smaller than the level spacing 5q even for T of the 
order of transition temperature T^.. 

The Hartree-Fock diagram as well as other diagrams, corresponding to single- 
particle level shifts and mixing, may be included as self-energy insertions into 
single-particle Green's functions in all diagrams of the previous subsections. 
As a result, the Green's functions become non-diagonal both in the orbital 
indices Z, and the spatial index p, if tunneling is involved. Such insertions al- 
ways introduce an additional smallness, as discussed in Sec. 6.5 for the case of 
tunneling in the ImSCBA. However, while inclusion of tunneling is necessary 
to go beyond the finite state space of a single localization volume, insertions 
non-diagonal only in the orbital index do not change the number of the final 
states and thus can be ignored. Insertion of tunneling into the renormalized 
interaction vertex, see Fig. 22a, produces interaction, nonlocal in space. How- 
ever, the corresponding correction to F/ is "small" compared to that of the 
ImSCBA diagram with Fig. 22b with the same final states, see Fig. 22c, d. Thus 
we conclude the variations of the Hartree-Fock potential do not generate any 
corrections to the transition temperature or the statistics of F; calculated in 
ImSCBA. 

However, there is an important feature in the statistics of the level shifts 
which changes the numerical factor in Eq. (140) and changes the power-law 
decay in the tail of the distribution of F;. To see this, let us consider now the 
second-order contribution shown on Fig. 4b. 2. Real part of this self-energy 
5^i{e) is included in SCBA but neglected in ImSCBA. To investigate its effect, 
we apply the Kramers- Kronig relation to F; '"■*(£) from Eq. (69a): 

5ii[e) = X5^ 2^ deide2de3F;^,i,.;3(ei,e2;e3)P— — -— ; 

h,l2,h ■' 6-61-62+63 

(156) 

where all the notation is introduced in Eq. (69a). Substituting Ai.{ei) = S{ei — 



77 




(e) 



~ \q\X^PT\n\\q\6cX^P) 



-iq 




~ \q\5(^\P\n{\q\5(^\P) 



-iq 

Fig. 22. Reduction of the non-SCBA diagrams (a-b) to the interaction vertex 
non-local in space (c). Characteristic function of the non-local interaction (c) in 
comparison with the corresponding block (f) of the SCBA diagram (d). 

^j.), we obtain 



hMz ^ ^'i + Sis 



whose characteristic function is 



|g|A(5c(T/T,). 



(157) 



(158) 



Therefore, the typical value of ~ \8q ^ b^. One thus might think that the 
main effect of this contribution is a weak random shift of the level position, 
so it can be disregarded as well. 

Nevertheless, despite its smallness, 5^/(e) introduces a new qualitative effect: 
repulsion between many-body levels. When S^i{e) is taken into account, one has 
to find the position of the shifted single-particle level from the equation 

e-6-<^6(e) = 0. (159) 

This equation describes repulsion between the single-particle excitation I and 
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three-particle excitations {h,l2',h)- solutions of Eq. (159) cannot approach 
each other by a distance smaller than A5^. 

In the considerations of Sec. 6 we assumed the energies of one-particle states 
to be independent of the energies of three-particle states into which these one- 
particle states decayed. As a result, resonant denominators [see, e.g., Eq. (123)] 
could become small independently of each other. The level repulsion suppresses 
such effect, so that Eqs. (130a), (132), and (135) overestimate the strength of 
the tail of the distribution. 

It is the same effect that was first discussed by Anderson [1] and analyzed 
rigorously by Abou-Chacra et al. [17]. According to their results, this effect 
leads to a change in the numerical coefficient in the equation for the critical 
disorder strength: e — > 2. Simply adapting this prescription for our Eq. (140) 
we obtain the transition temperature modified by level repulsion 

12XMTc , 1 

In - = exp 

A 

Now we will sketch the derivation of Eq. (160) for our problem. First, we have 
to identify the sequence of the diagrams which may give the level repulsion 
between the resonant multiparticle states. To do so, we include the level shift 
S^i{e) into the imaginary part of the fourth order SCBA diagram as shown on 
Fig. 23a. 

Direct inspection of the diagrams Fig. 23a shows that for the levels , . . . , 
maximizing the skeleton diagram (aO) of Fig. 23, only diagrams (al) and 
(a2) give rise to the simultaneous divergence of 6^, whereas in the remaining 
diagrams (a3)-(a6) 6^ fluctuates independently of F^^. 

It means that the role of the level shift 5^ in the diagrams (a3)-(a6), is, indeed, 
just the broadening of the distribution of C,i by the value of \5q{T /T^) <^ 5^, see 
Eq. (158). This broadening does not significantly change the distribution of the 
resonant denominators and can be disregarded together with the fluctuations 
of the Hartree-Fock potential. On the other hand, the level shifts on diagrams 
(al), (a2) are large whenever F^^ is maximal and they describe the effect of 
level repulsion discussed before. 

The maximally divergent series incorporating the correlations in and F; 
is generated by the hnearized self-consitent-Born approximation (LSCBA) 
shown on Fig. 23b, as it can be checked by explicit consideration of several 
iterations. 

The analytic expression of the LSCBA self-energy is [c/. Eqs. (69a), (105a))]: 



Ci5q |ln /| 



(160) 
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rK6) = rr^H6)+i: 



P5pA[e-^iAP + 0')]^i{e:P + a) 



t% [e-6i(p + a) -(^6i(e;P + a)]^ 



h,i2,i3 



+ 



9/s,{e - E\li^)riAEi^i^ - e) 



> . 



(161b) 



As in Eq. (69a), the coordinate p is assumed to be the same in all terms 
in the equations unless it is specified explicitly otherwise. The notation HJ^^ " 
was introduced in Eq. (113). The ultraviolet cutoff function 9a{x) is the same 
as in Eq. (104). It is important to emphasize that it depends on the unshifted 
energies of the levels. This is because the cutoff was introduced in the non- 
interacting Green function first; the self-energy appears in the denominator of 
the SCBA-dressed Green function as a result of the summation of a geomet- 
ric series; such summation changes the denominator only, keeping the cutoff 
intact. 

Now let us go through the steps of Sec. 6.4 and see how they are affected by 
the shifts in the denominators of Eqs. (161). First, we notice that Eq. (161b) 
is still given by the diagrammatic representation of Fig. 11 with the change in 
the rules of reading the double fine: 

e ^_^A(e-6) 



[e-Ci-SCiie] 



l2 • 



(162) 



The lowest order term of the cluster expansion, Eq. (114), see also Fig. 13, 
becomes 



p(i,o) ^ 2r(°'°) 



_^p(0,0) j2 



1,^3 [e-E\li^-d^i,{e-Ell 



(163) 



/i,i2,;3 [e - E\li^ + 5^/3(2^1^2 - e) 
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We see that the corresponding potential is no longer three-particle, but de- 
pends on all the coordinates via 6C,, which depend on the positions of all 
other levels. However, when performing the linked cluster expansion, in each 
term of the sum over triples {hjhjs) we can shift an integration variable: 
61 + ~^ 61 ■ After this shift Eq. (15) acquires the form 



"■°'=in(/f E 



/l2: 



;3 



^ 1=1 



( „;=ii2cosh|, 
1) {0^[e-El, + 5^,)y, 



ir+n, 



(164) 



A3 = (e-^^^-^ -1) (^6-S?,-56))', 



where (...)' denotes the average (73) with levels 1,2,3 excluded, and the 
potentials U^l'.^^ were defined in Eq. (115). Using the fact that 6^i{e — Sj^jj) 
from Eq. (161a) is not singular at — > e and \5^\ <^ A, we can perform the 
integration in the same way as in Sec 6.4 to obtain 



(165) 



Here the only effect of the shifts was to make the electron contribution statisti- 
cally independent from the hole one. As a result, the cross-term in Eq. (116b) 
vanishes and the numerical coefficient is changed in comparison with Eq. (117). 

Let us now consider the modification of the fourth order result (123). After 
two iterations of Eq. (161b), see also Fig. 14 and Eq. (162), we find 



X 



+ 
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+ 



hk 
2UI5 



_ _ 



e — ^ 



me 



e - - 56,(e - E't,)] [e - Sj^je^^ + 5^1^ (e\1,^,^ - e) 



(166) 



where the prime in the second sum has the same meaning as in Eq. (123). Each 
term in Eq. (166) produces its own cluster function, see Eq. (125). Similarly 
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Fig. 23. (a) Including the level shift into the fourth order ImSCBA diagram. Only 
diagrams a.l and a. 2 are resonant, i.e., divergences in 6(,i and Ti are correlated, (b) 
Linearized SCBA approximation which incorporates the simultaneous divergences 
in S^i and F;. 

to Eq. (164), the cross-correlation terms vanish and we obtain after obvious 
shifts of the variables 
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5ei(Si2 - e) 
<^6(5i2 - e) 



^^A 



e — 



77126 



6 - Etl, + SCe{e - Etl) 



(167c) 

where (...)' denotes averaging (73) with the levels 1, 2, 3, 4, 5, 6 excluded. The 
potentials U here are defined in Eq. (124). 
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Apparently, the role of the level shifts in Eqs. (167c) is similar to that in 
Eq. (164) - perturbative modification of the cut-off. However, unlike the lowest 
order perturbative corrections, the shifts in Eq. (167c) contain a resonant term 
which depends on variables Ci, ■ ■ ■ ,^6 only and therefore can not be treated 
as a non-correlated random number. To see the origin of such resonant term, 
consider the argument of the first ^-function in f^"^. Prom Eq. (161a) with 
/ = 0, we find after the same variable shifts as in /^^ 



S^i{e - El) = 



^ ~ --452 



+ ..., 



(168) 



where . . . denote the terms which contain extra levels. Those terms are random 
and can be disregarded. One can check by the same method that the shift 
5^^€ — SI4) always depends on extra levels and does not produce resonance. 

Having in mind that the ^^-functions in Eqs. (167c) cut off the logartihmic 
divergence, we neglect the non-resonant terms, and simplify Eqs. (167c) as 
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(169) 



Substituting Eq. (169) into Eqs. (167a)-(167b) and performing integration we 
obtain analogously to Eq. (127) 

lnW^(2,o)^_(^^^r(o,o)^i 

In addition to the change in the overall numerical factor similar to Eq. (165), 
Eq. (170) shows new important feature. Namely, the argument of the logartihm 

is no longer dependent on the parameter s, or, in other words, the algebraic 
tail in the distribution function is suppressed in comparison with Eq. (127). 
This is the manifestation of the level repulsion in complete analogy with the 
arguments of Anderson [1] . 

^® Keeping them will be beyond the accuracy of LSCBA approximations where the 
terms of the same order were neglected in the very beginning. 
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Fig. 24. The cluster function j^^^ for the six order perturbation theory expansion. 
The notation is defined on Figs. 11, 12 and in Eq. (162). 



The procedure outlined above can be continued to all orders of perturbation 
theory. In particular, for the sixth order cluster function. Fig. 24, one finds 
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(171) 



which means that the scale of the integration is determined by the energy 
of the previous generation only. This transfer matrix structure of the cluster 
functions repeats itself in all orders and makes it possible to perform the 
integration in any order. Instead of Eq. (135) we find 
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C(n, m) 



/In 



(172) 



i.e., once again the algebraic tail becomes "more short-range" due to the 
suppression of the large denominators. Coefficients C{n,m) are insensitive to 
the resonant denominators and they are still given by Eq. (136). We, therefore, 
can repeat all of the arguments after Eq. (136) and obtain Eq. (160) after 
replacement (128). 



To conclude, the level repulsion present in SCBA but neglected in SCBA leads 

to the change in the numerical prefactor in the expression for the transition 
temperature [cf. Eq. (160) with Eq. (140)] but does not affect the statement 
about the existence of the transition itself. All other corrections lead to the 
perturbative corrections to T^. 
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7.^ Validity of ImSCBA in the metallic phase 

In the metallic phase the classification of non-SCBA diagrams into interaction 
vertex corrections and interference terms is the same as presented in Sees. 7.1 
and 7.2 for the insulating phase. What changes is that the resonant terms in 
all of the expressions now acquire the finite width F*^*"^) 

r r*^*") 1 X 

•KSix) + ^; P > ^, (173) 

where F^*'*-' is given by Eq. (74). This finite width prohibits vanishing of the 
denominators, thereby cutting off the power-law tails of the distribution func- 
tions for all the quantities considered in Sees. 7.1-7.3. As a result, all the 
cumulants of the distribution function become finite. Moreover, under the 
condition (80) the distribution functions may be considered Gaussian. 

Another feature is that the level occupation (50) is no longer binary, ni{e) — 
±1, as in insulating phase, but it is kept close to its equlibrium value ni{e) — 
tanh by the energy relaxation. This further suppresses the fluctuations by 
the factor of 5c^/T. 

On the other hand, in Sec. 7.1, we established that the role of the higher-order 
correction to the vertices is the perturbative renormalization of the SCBA 
results. As the scale of the fluctuation of those vertices is suppressed even 
further in the metallic phase, neglecting this renormalization is justified even 
more in this phase. 

The role of the particle permutations and interference in the final state, see 
Sec. 7.2 is investigated by direct evaluation of the mesoscopic fluctuations and 
averages of the diagrams of Fig. 19 in a fashion of Fig. 8. Those contribution 
are smaller than SCBA values by a factor of M at least. Therefore, deep in 
the metaUic phase the kinetic equation consideration of transport in Sec. 5.2 
is well justified. 

Finally, in the insulating phase the most important effect, not included in 
ImSCBA, was many-mody level repulsion, see discussion around Eq. (159). 
The energy scale of this repulsion is \5q. On the other hand, at T > Tj„, where 
Tin is defined in Eq. (80), we find F^^") > \5^. This means that the spectrum 
structure produced by this repulsion is smeared by inelastic processes and 
need not be taken into account. 

Thus, we conclude that the ImSCBA is justified parametrically in deep in the 
metallic phase. 



85 



8 Conclusions 



This paper is devoted to the analysis of the low-temperature transport in 
disordered conductors with localized one-particle states. The main question 
is whether the electron-electron interaction alone is sufficient to establish the 
thermal equilibrium in the system. The same question can be formulated even 
more boldly - whether there is a many-body mobility threshold, i.e., energy 
separating the many-body states localized in the Fock space of the system 
from the states which are delocahzed. 

One can apply weak-localization arguments to show that conductivity at high 
enough temperatures is non-zero. It is not disputable as long as the inelastic 
dephasing rate l/r^ exceeds the level spacing in one-particle localization vol- 
ume (cell) 5(^. Extension of this approach to lower temperatures is problematic 
as the quantum corrections to conductivity diverge. For this reason, in order 
to describe the low-temperature behavior, we adopt a different strategy. We 
formulate the effective Hamiltonian description for the processes with the en- 
ergy transfer of the order of 5^. Reduction of the original Hamiltonian to the 
effective one is not performed systematically. Nevertheless, we beheve that it is 
an apropriate low-energy limit of the theory of electrons in disorder potential. 

Statisitical analysis based on the effective Hamiltonian enables us to demon- 
strate the stability of two qualitatively different phases - metallic, for T > T^. 
and insulating for T < T^, where Tc is given by Eq. (160). This corresponds 
to the existence of the many-body mobility threshold £c related to Tc by the 
thermodynamic formula (101). 

We show that deep in the metalUc phase, T ^ (see Sec. 5) the trans- 
port coefficients in the system are self-averaging. Using this fact, we derived 
the quantum Boltzman equation. The temperature dependence of the elec- 
trical conductivity o'(T), following from this equation (see Sec. 5.2), is quite 
non-trivial even for T ^ Tc. Namely, (t{T) increases with T as at low 
temperatures, while at high temperatures it saturates, approaching the Drude 
limit. Thermal conductivity deviates from the Wiedemann- Franz law with the 
decreasing temperature. However, this deviation is never strong, see Fig. 10. 

In the insulating phase, T < Tc, we use Feynmann diagram technique to 
determine the probability distribution function for quantum-mechanical tran- 
sition rates. The probability of an escape rate from a given quantum state 
to be finite turns out to vanish in every order of the perturbation theory in 
electron-electron interaction. Thus, in the absence of coupling to any external 
bath (phonons) electron-electron interaction alone is unable to cause the re- 
laxation and establish the thermal equilibrium. In other words, the insulating 
phase is stable. 
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Although a{T) — exactly as long as T < Tc, the stability of the insulator de- 
creases as T approaches T^. It means that effects of interaction of the electrons 
with the external bath (phonons) become more and more pronounced. More 
precisely, if the electron-phonon coupling is weak and T <^ Tc, one phonon 
can cause at most one electron to hop because the phase volume of the acces- 
sible final states is quite small. The closer is T to Tc, the bigger is this phase 
volume, and in the vicinity of the transition point one phonon can initiate 
a whole cascade of the electronic hops. The size of the cascades fluctuates 
strongly, depending on the realization of disorder in the system. The typical 
size of the cascade grows at T —> Tc [see Eq. (142)]. It means, that even the 
inflnitesimal electron-phonon interaction would produce a finite conductivity. 
This is the onset of the metallic phase. 

We also conclude that the phonon-induced hopping conductivity in the in- 
sulating phase close to the transition is strongly enhanced by the electron- 
electron interaction. This conclusion can be relevant for the numerous experi- 
ments [32-36] , where the observed conductivity in the strongly localized phase 
of disordered conductors was too large to be explained by conventional theory 
of phonon- assisted hopping conductivity. 

It should be emphasized that the many-body localization, which we discuss 
in this paper, is qualitatively different from conventional finite temperature 
Metal to Insulator transitions, such as formation of a band insulator due to the 
structural phase transition or Mott-Hubbard transition [37,38]. In these two 
cases, at a certain temperature T* a gap appears in the spectrum of charge 
excitation (Mott insulator) or all excitations (band insulator). However the 
conductivity remains finite although exponentially small as long as T > 0. 
This is not the case for many-body localization, which causes exactly zero 
conductivity in the low-temperature phase. 

Is the many-body localization a true thermodynamical phase transition with 
corresponding singularities in all equilibrium properties? This question defi- 
nitely requires additional studies, however, some speculations can be put for- 
ward. The physics described in the present paper is associated with the change 
of the characteristics of the many-body wavef unctions. It is well known that 
for non-interacting systems localization-delocalization transition does not in- 
fluence the average density of states, i.e., it does not affect any macroscopic 
thermodynamic properties. Application of the same logic to the exact many- 
body eigenvalues would indicate that the many-body localization tranisition is 
not followed by any singularities in the static specific heat, etc. On the other 
hand, at this point we can not rule out the possibility that this conclusion 
is an artefact of treating the real parts of the electron self-energies with an 
insufficient accuracy. Most likely scenario, to our opinion, is that the insulat- 
ing phase behaves like a glass (spin or structural) and demonstrates all the 
glassy properties [19], like absence of ergodicity (even when some coupling 
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with phonons is included), effects of aging, etc. Question of the equihbrium 
susceptibihties in the latter case becomes quite meaningless. 

The quantitative theory built in this paper assumes that the interaction is 
weak. On the other hand, qualitative consideration of the localization of many- 
body excitations does not rely on this assumption. The important ingredients 
are (i) localization of single-particle excitations, and (ii) Fermi statistics. Con- 
sider, as an example, Wigner crystal [39]. It is well known that strong enough 
interaction leads to a spontaneous breaking of the translational symmetry in 
d-dimensional clean systems at d > 2. In a clean system Wigner crystaUiza- 
tion is either a first-order phase transition (rf = 3), or a Kosterlitz-Thouless 
transition {d = 2). Even weak disorder destroys both translational and orien- 
tational order [40] and pins the crystal. The symmetry of this state is thus 
not different from the symmetry of a liquid, and the thermodynamic phase 
transition is commonly believed to be reduced to a crossover. 

We argue that the many-body localization provides the correct scenario for the 
finite-temperature "melting" transition between the insulating phase, which 
may be called "solid" , and the metallic phase, which may be called "liquid" . 
Indeed, the conductivity of the pinned Wigner crystal is provided by the mo- 
tion of defects. At low temperatures and in the absence of the external bath, all 
defects are localized by the one-particle Anderson mechanism. Phonon modes 
of the Wigner crystal are localized as well, so the system should behave as a 
many-body insulator. As the temperature is increased, the many-body metal- 
insulator transition occurs, though it is not clear at present, whether it occurs 
before or after the crystalline order is destroyed at distances smaller than 
Larkin's scale. Construction of effective theory of such a transition is a prob- 
lem which deserves further investigation. 



9 Acknowledgements 

We are grateful to M. E. Gershenson and V. E. Kravtsov for valuable remarks. 
We also want to thank A. D. Mirlin for clarifying statements made in their 
preprint [28] and reveahng the content of Eqs. (19)-(20) of Ref. [28] to us. 

I.L.A. was supported by Packard foundation, B.L.A. acknowledges financial 
support from US DOE Office of Science under contract No. W-31-109-ENG- 
38, ARO/ARDA (DAAD19-02-1-0039) DARPA under QuIST program and 
NSF grant DMR 0210575. 

Note added- When the current manuscript was almost completed, a new ver- 
sion of preprint [28] appeared [41], which deals with some of the problems 
constituting the content of the present paper. 



88 



A Slow energy relcixation? 



Let us now come back to the assumption of the strong inelastic relaxation 
which was used to validate expansion (91). One could think that in some 
cases this assumption is not valid, which would lead to the deviation of the 
temperature dependences from those of Eqs. (93). To investigate the limits 
of validity of Eqs. (93)-(100), we employ the following qualitative arguments. 
These arguments only slightly modify the discussion after Eq. (99). 

We notice that the transport still occurs through rare pin-holes. However, the 
insufficient rate of the inelastic processes allows the electron to explore an 
energy strip of the width e <^ T before leaving the site via tunneling to the 
neighboring sites. Correspondingly, instead of using Eq. (91) one may look for 
the solution in the form 

$(6, p; t) = P{e/e) , / dxPix) ^ 1 , (A.l) 

where e is a scale to be found self-consistently. The conductivities are still 
given by Eqs. (93), but the random quantities 8^,3^ replaced by different 
ones BctjBk, which are defined through f3{xT/e) rather than through /9o-,«;(a:). 
Thus, either a straightforward calculation, or a qualitative argument, similar 
to that after Eq. (99), give for the typical B an estimate 



p(in) 



e 



B ^ , (A.2) 

valid \IB 1. To find e, we balance the elastic and inelastic terms in Eq. (90). 
The estimate for the elastic term is 

PS"^ ^(^(e, p) A{e, p + a))$ ~ Pd^m ~ T^^") — $, (A.3) 
a (><: 

while the inelastic term is estimated from Eq. (87c) as: 

Requiring these two rates to be of the same order, we estimate the scale e as 

6- = 5, (- j . (A.5) 



Conditions e <^ T and B <^ 1 would produce a parametric temperature region 
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in which the non-equihbrium function <l>(e) does not have the quasi-equihbrium 
shape (91). However, existence of such a regime requires the condition 



//\4/3 / 1 \7/6 , 



(A.7) 



which is not consistent with the additional assumptions M^A ~ 1, / ~ A <^ 1. 
Therefore, the regime of the slow energy relaxation is not feasible for this model 
and Eqs. (93), (94) describe the entire temperature dependence for T > Tin- 



B Probability distributions for Ba,K- 



Let us represent the definitions (92c) in the form 



)2 + ' 



(B.l) 



where the functions f3a{x) and Pk{x) are given by Eq. (92d). (We omit the 
superscript "(in)" of F*^™) everywhere in this appendix). The positions of the 
levels are assumed to be completely uncorrelated, so the sought charac- 



teristic function P(s) 
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can be represented as 



P(s) — lim 
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(B.2) 



Equation is equivalent to the partition function of the classical gas of two 
species, (±), and interacting to each other via pair potential B{^^,^^). It 
can be immediately evaluated using Mayer-Mayer cluster expansion [29]. We 
will keep the contributions upto four-particle clusters to justify the further 
approximations : 
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where 



/i2 = exp[-s5(6,6)]-l, (B.3e) 

is the Mayer's /-function which vanishes rapidly with the distance between 
the levels |^i — G|. Performing integrations in Eqs. (B.3b)-(B.3d) with the 
help of definition (B.l) and using F ^ T, we find 
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(B.4) 



Here Xq and Xi are the modified Bessel functions. We also have introduced the 
parameter r = SttFT/S'^, see Eq. (96). 

Next step is to notice that T ^ 5^ and the observable quantities are con- 
tributed by s such as rS2{s/r) < 1. Therefore, the contribution from the 
higher cluster are suppressed and can be neglected. All the further calculation 
is performed using P^^-* only. 



Two hmiting cases can be considered. First, if r » 1, S{x) can be expanded 
S{x) — x{l — x/2 -I- . . . ), which gives 



r > 1 : <^ 
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(B.5) 



In the opposite case, r <^ 1, we use the asymptotic expansions of the Bessel 
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functions to approximate S{x) ~ ^J2/{t:x), to obtain 
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where G ^ 0.916 is the Catalan's constant. 

It is instructive to calculate the inverse Laplace transform of P^i^[s) and find 
the distributions functions Po-,k('B). In both limiting cases of large and small r 
one can use the saddle point approximation. For r 3> 1 the distributions are 
approximately Gaussian with {Ba,K) = 1 and 
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corresponding to self-averaging of Ba,K- For r <S 1 have: 
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It is noteworthy that, for B ^ r this ditribution coincides with the distribution 
of the largest of n = VSttc {T/5c) independent random variables 



26iT/-K 
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4T a '■■■'4r a 

[see Eq. (92c)], with uniformly distributed in the range < < 5^/2. 



(B.9) 



C Cancellation of the cutoff A 



The formal reason to introduce the cutoff A in Sec. 6 was the insufficient ability 
of the 3n-particle potentials to confine particles, which resulted in logarithmic 
divergences in the leading terms of the cluster expansion for n > 1. On the 
other hand, since the transition is associated with anomalously small values 
of the energy denominators, rather than anomalously large ones, it is natural 
to expect the transition condition not to contain A at all. This means that A 
should cancel out when higher-order terms of the cluster expansion are taken 
into account. Moreover, this cancellation should occur for each component of 
the potential separately {e.g., for each of the 6-particle potentials C/^^, C/^'*, 
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Fig. C.l. The diagrams canceling the cutoff A in the fourth diagram of Fig. 14. 



C/'*^, and U'^^). In this appendix we first show how the cancellation occurs for 
the 6-particle potential U^^, and then discuss 3n-particle potential JJ^^-^ for 
an arbitrary n. The terms corresponding to an arbitrary sequence of electron 
and hole decays can be analyzed analogously. 



Consider connected terms of the cluster expansion containing k functions 
f^^{^i\ ■ ■ ■ , d^^) V ■ ■ , f^^{^i\ Ce^)- To maximize the corresponding con- 
tribution one should connect them in such a way that one of the energy de- 
nominators coincides for all terms (two denominators cannot coincide, as it 
would require the coincidence of ^P'* for all j and each i, which is prohibited 
by the rules of the cluster expansion). Having fixed one of the denominators, 
one should keep the remaining free level positions independent of each other in 
order not to obtain smallness in M or T/S/^. As a result, two diagrams should 
be evaluated for each k, shown in Fig. C.l. Summing over k and adding the 
leading term, we obtain 



k k\{k-i) , 



(2XMT\' 



In 



dx cosh^^^ ~ f 2XMT \ 



k-l 



COS! 



(C.l) 



where the logarithmic term in the square brackets is the leading one, corre- 
sponding to the fourth diagram of Fig. 14, the second and third terms corre- 
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spond to the diagrams of Fig. C.l. The summation over k is performed as 

oo 



(C.2) 



l/A 

As a result, we obtain 



This seems to correspond to the replacement 



A^-yMT^, (C.4) 
rather than to the promised Eq. (128). 

To show that the transition point is nevertheless determined by Eq. (128), 
we analyze the hole channel for an arbitrary n. In this case the sublcading 
logarithmic correction is given by diagrams where n — 1 energy denominators 
are fixed by connections, and only one is independent of others. The indepen- 
dent denominator can be chosen in n ways, of which n — 1 produce diagrams 
analogous to the first diagram of Fig. C.l (each / has only one argument inde- 
pendent of other /'s) and the last denominator produces a diagram similar to 
the second diagram of Fig. C.l (each / has three last arguments independent 
of other /'s). The result is 



X 



(n-1)! ^/s^(M (n-2)! VsfW V k 



(C.5) 



The expression in the square brackets is nothing else but the first two terms 
of the binomial 



1 r (A/XScT , /MTA"^'" ^ 
m — , — m ' 



(n-1)! [ ^/s^m \ 5^ 52 



(C.6) 



(n-1)! V^^A" Vsr(O'O), 



To obtain the next terms of the binomial one would have to connect the cluster 
functions in a way that fixes n — 2, n — 3, etc. denominators. As seen from 
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I 



Fig. D.l. Amplitude of nth order of the peturbation theory used for the estimate 
of the distribution function (D.l). Being squared, these amplitudes generate the 
electron-like decay processes only (the first diagrams on Figs. 13, 14). External lines 
are assumed to be amputated and shown for the notation of orbitals and energy. 
Path £ in Eq. (149) corresponds to the sequence {^1,^2, •• • ^n-i} on this figure. 

Eq. (C.6), the factor under the logarithm becomes unimportant at 

large n. This justifies the rule (128). 



D Derivation of Eq. (151) 

Our goal in this Appendix is to evaluate the most dangerous contribtion from 
the interference corrections to the SCBA result. For this purpose, it suffices 
to demonstrate the calculation of from Eq. (151) without numerical coe- 
ficient. Therefore, we will take into account only electron- like processes shown 
on Fig. D.l and disregard other processes, cf. Figs. 13, 14. This leads to the 
underestimate of the overall numerical coefficient but does not affect the fac- 
torials. 

The expression for such nth order amplitude in the insulating region reads 



A^({4};{/fc};W}) 




In-" f,l 



^1 jl2---iln i^n—l 



(D.l) 



where the notation was introduced in Eqs. (105b) and (113). 
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Expression (D.l) has the same form as those considered in Sees. 6.4 and 7.1. so 
that the machinery of Sec. 6.3 can be apphed with small modification. First we 
fix the electron and hole levels C/d • • • > C/„j imn ■ ■ ■ , Cm„ in Eq. (D.l), and sum 
over all the intermediate states £i, £2; • • • , ^n- Calculating the characteristic 
function, we perform the partial averaging - i.e., procedure of Eq. (102) with 
2n levels Zi, . . . , Z„; mi, ... , m„ excluded. 

Using the linked cluster expansion and Eq. (48g), we find 



In W {q; {k}; {rrii}) = In ^exp 



i{^k};{ikh{mk}) 



{fc,};{mfe} 



-7r|g|(5^A 



2Alni 
A 



k=l 



rrik 



(D.2) 



r^imim2...mn-i 

In 



where we were dealing with the logarithmic integrals as in Sec. 7.3. For the 
sake of brievity we replaced S;"^" /5^ — >• in the last two lines of the 

equation. 

Equation (D.2) enables us to estimate the distribution function of the am- 
plitudes summed over the permutations among n electrons and n holes, see 
Eq. (149). Approximating the amplitudes obtained by the permutations of the 
final states to be independent of each other, cf. Fig. 20 we find 



W{q;{h};{mi}) = ( exp 



iq {-lf'+''-A{{Pi{k};Pm{m,}) 



PlPm 



(D.3) 



n W{q;Pi{k};Prr^{mi}), 

PlPm 



where we used the standard expression for the characteristic function of sum, 
we suppressed the ik arguments of the amplitude for brevity, and the partial 
averaging as in Eq. (D.2) is imphed. 



This corresponds to the leading logarithmic approximation in the Unked cluster 
expansion. Note, however, that it does not affect the terms divergent as n! at all. 



96 




7/({4},{^fe}) 

{lk},{mk} 

.e-^TT— V P ^ (D-4) 
W(g;{l,2,...n};{n + I,...,2n})-1 



X 



(n!)^ 



Finally, expanding the last in Eq. (D.4) in powers of |g|, and performing re- 
maining integrations, we obtain the structure of Eq. (151). Indeed, exponent 
contains (n!)^ terms itself, so that the factorial factors cancel at all in the lead- 
ing term. In the next term product W (q'; {Z^}; {mfe}) W (5; Pi{Zfc}; P^lmfc}) 
gives the same volume for the ^ integration only for Pi = and impose the 
additional restricitions otherwise. The factor of n! is, thus, just a number of 
permutations electron-hole pairs without their destruction. 
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